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Abstract 



We present the implementation of the radiative corrections of the Higgs sector in three 
public computer codes for the evaluation of the particle spectrum in the Minimal 
Supersymmetric Standard Model, SoftSusy, SPheno and SuSpect. We incorporate 

i-S3 the full one-loop corrections to the Higgs boson masses and the electroweak symmetry 

breaking conditions, as well as the two-loop corrections controlled by the strong gauge 
coupling and the Yukawa couplings of the third generation fermions. We include also 

_C ■ the corrections controlled by the r Yukawa coupling that we derived for completeness. 

The computation is consistently performed in the DR renormalisation scheme. In a 
selected number of MSSM scenarios, we study the effect of these corrections and analyse 
the impact of some higher order effects. By considering the renormalisation scheme and 
scale dependence, and the effect of the approximation of zero external momentum in 
the two-loop corrections, we estimate the theoretical uncertainty on the lighter Higgs 
boson mass to be 3 to 5 GeV. The uncertainty on M^ due to the experimental error 
in the measurement of the SM input parameters is approximately of the same size. 
Finally, we discuss the phenomenological consequences, using the latest value of the 
top quark mass. We find, in particular, that the most conservative upper bound on 
the lighter Higgs boson mass in the general MSSM is M^ < 152 GeV and that there 
is no lower bound on the parameter tan (3 from non-observation of the MSSM Higgs 
bosons at LEP2. 
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1 Introduction 

The Minimal Supersymmetric Standard Model (MSSM) pQ provides an attractive weak-scale 
extension of the Standard Model (SM). At the price of introducing a bosonic (fermionic) 
super-partner for every fermionic (bosonic) particle of the SM, plus an additional doublet of 
Higgs fields and their fermionic partners, it allows one to solve the hierarchy problem, leads 
to a consistent unification of the SM gauge couplings and provides a natural candidate for 
the dark matter in the universe. In the last two decades, a large effort has been devoted 
to the detailed theoretical study of the properties of the supersymmetric particles and the 
MSSM Higgs bosons, and to the experimental search for these particles, either directly in 
high-energy collider experiments or indirectly through high-precision measurements. 

The task of providing a detailed phenomenological analysis of the masses and couplings of 
the supersymmetric particles and the MSSM Higgs bosons, and comparing it with the results 
or expectations of the present and future experiments, proves to be an extremely difficult one, 
due to the large number of new parameters (>100 in the most general soft supersymmetry 
(SUSY) breaking Lagrangian and >20 when some phenomenological constraints are imposed; 
see Ref. |2] for a discussion). However, there are well motivated theoretical models where the 
soft SUSY-breaking parameters obey a number of universal boundary conditions at some 
very high scale, such as the Grand Unification Theory (GUT) scale Mgut, or the messenger 
scale M mess , thus reducing the number of fundamental parameters to a handful. This is 
for instance the case of the models of gravity mediated (mSUGRA) [3], gauge mediated 
(GMSB) [4J or anomaly mediated (AMSB) |Hj SUSY breaking. String-inspired models of 
soft SUSY-breaking have also been proposed 0. In all such models, the soft SUSY-breaking 
parameters are given as input at the high-energy scale, then they are evolved by means of the 
MSSM Renormalisation Group Equations (RGE) down to a low-energy scale, Mewsb, where 
the conditions of electroweak symmetry breaking (EWSB) are imposed and the spectrum of 
superparticle masses and couplings is computed. 

Several computer codes (that in what follows we will often denote as spectrum generators) 
have been developed in the past years to provide reliable determinations of the supersymmet- 
ric spectra in models with high-energy boundary conditions, including the calculations of the 
various radiative corrections. In particular, the publicly available codes include IsaJet jjj, 
SoftSusy jH], SuSpect jH] and SPheno [TU] . Detailed comparisons of the features and the 
results of these codes can be found in Ref. ^1] for instance. 

A crucial prediction of the MSSM is the existence of at least one light Higgs boson ^2] 
which, at the tree level, is bound to be lighter than the Z boson. If this upper bound was not 
significantly raised by radiative corrections, the failure of detecting this Higgs boson at LEP 
would have ruled out the MSSM as a viable theory for physics at the weak scale. However, it 
was first realised in Ref. J3] that the inclusion of the one-loop 0(a t ) corrections, which rise 
quartically with the mass of the top quark and logarithmically with the mass of its scalar 
superpartner, may push the lighter Higgs boson mass well above the tree-level bound. 



In the subsequent years, an impressive theoretical effort has been devoted to the precise 
determination of the Higgs boson masses in the MSSM. A first step was to provide the full 
one-loop computation, performed in Refs. fU E3 E] ■ A second step was the addition of 
the dominant two-loop corrections which involve the strongest couplings of the theory: the 
QCD coupling constant and the Yukawa couplings of the heavy third generation fermions 
(although the masses of the bottom quark and the r lepton are relatively tiny compared to 
the top quark mass, the b and r Yukawa couplings can be strongly enhanced for large values 
of tan/?, the ratio of the vacuum expectation values of the two Higgs doublet fields). 

The leading logarithmic effects at two loops have been subsequently included via appro- 
priate RGEs [TH EU, and the genuine two-loop corrections of 0(a t a s ) [T3 ] I2TI ] I2T ] 122 ] 123] 
and 0(ab<y s ) [23] have been evaluated in the limit of zero external momentum. The two-loop 
Yukawa corrections of 0(af) [TRESES] and 0(a t ab+al ) [2E] have been also evaluated in the 
limit of zero external momentum. To complete the calculation of the two-loop corrections 
controlled by third-generation fermion couplings, only the expectedly small corrections that 
are proportional to the r-lepton Yukawa coupling are still to be determined. The 0(a^) 
corrections can be straightforwardly derived from the 0(af) corrections as explained in 
Ref. |25] ; the missing C(a fe a T ) contributions are provided in this paper. 

The tadpole corrections needed to minimise the effective scalar potential, V e {[, have also 
been calculated at the one-loop (TBI EZ| and two-loop [2H1 EHj levels for the strong cou- 
pling and the top and bottom quark Yukawa couplings; the tadpole corrections involving 
the r-lepton Yukawa coupling will also be provided in this paper. Finally, the full two- 
loop corrections to the MSSM effective potential have been calculated [2H], together with 
a first study of the two-loop corrections to the lighter Higgs boson mass controlled by the 
electroweak gauge couplings (30|. 

The calculation of the radiative corrections to physical observables requires the choice of 
a renormalisation scheme. For example, one might choose to express the corrections in terms 
of "On-Shell" (OS) parameters, such as pole particle masses and suitably defined mixing 
angles. This is the scheme adopted for the computation of the MSSM Higgs boson masses 
by the computer code FeynHiggs [2], based on the results of Refs. [THIEO] and subsequently 
expanded to include the results of Refs. [2111231201122] (see Ref. [3S] for a discussion). 

However, when the MSSM parameters at the weak scale are derived from a set of unified 
parameters at the GUT scale via renormalisation group (RG) evolution, they come naturally 
as unphysical "running" quantities expressed in a minimal subtraction scheme such as the 
modified Dimensional Reduction scheme, DR, which is usually adopted since it preserves 
SUSY, at least up to two-loop order. Interfacing the output of the RG evolution of the 
MSSM parameters with a two-loop calculation of the Higgs boson masses based on the OS 
scheme requires translating the DR input parameters into physical quantities, at the two- 
loop level for the parameters entering the tree-level mass matrix (such as the CP-odd Higgs 
boson mass) and at the one-loop level for the parameters entering the one-loop corrections 
(such as the squark masses and mixing angles). This induces additional complications in the 
procedure. 



A more direct strategy would be to perform the computation of the Higgs boson masses 
directly in the DR renormalisation scheme, using the MSSM parameters as they come from 
the RG evolution. The results for the Higgs boson masses must be equivalent to those 
of the OS calculation up to terms that are formally of higher order in the perturbative 
expansion. The numerical differences between the results of the two calculations can be 
taken as an estimate of the size of the corrections that are still uncomputed. They can be 
viewed, together with the differences due to changes in the renormalisation scale at which 
the corrections are evaluated, as part of the theoretical uncertainty in the calculation. 

A first purpose of this paper is to present the implementation of a purely DR calculation 
of the neutral MSSM Higgs boson masses, based on the one-loop results of Ref. fB] and the 
two-loop results of Refs. [231 123 EH1 EH] , into the latest versions of three public codes for the 
RG evolution of the MSSM parameters and the computation of the superparticle and Higgs 
boson mass spectrum, i.e. SoftSusyl .8.7, SPheno2.2.1 and SuSpect2.3. We also derive 
the small 0{a\ ) a T ) contributions to complete the calculation of the two-loop corrections to 
the Higgs boson masses controlled by the third-generation Yukawa couplings. 

A second purpose of this paper is to estimate the various theoretical uncertainties in the 
determination of the Higgs boson masses. Using representative choices for the high-energy 
boundary conditions, known as Snowmass Points and Slopes (SPS) [31], we first compare 
the outputs of the three codes, which are generally in very good agreement. We outline 
the reasons for the small residual discrepancies and briefly discuss the effect on the Higgs 
boson masses of various two-loop refinements of the general calculation. We then compare 
the results of the DR computations of the Higgs masses with those of the OS computation 
in the program FeynHiggs and obtain an estimate of the residual scheme dependence. We 
also discuss the renormalisation scale dependence of our results for the physical masses, 
which can also be taken as an estimate of the theoretical uncertainty, and the effect of the 
approximation of zero external momentum in the two-loop corrections. A similar analysis is 
performed in the context of the unconstrained MSSM. Finally, we estimate the uncertainties 
in the Higgs boson masses arising from the experimental error in the determination of the 
top and bottom quark masses and the strong and electromagnetic coupling constants. 

A final purpose of this paper is to study the impact of these results on MSSM Higgs 
phenomenology, under different assumptions for the mechanism of supersymmetry breaking. 
In particular, we estimate the upper bounds on the lighter Higgs boson mass in the general 
MSSM and in the mSUGRA, GMSB and AMSB scenarios once the latest values of the 
SM parameters (in particular the new experimental value for the top quark mass, M t = 
178.0±4.3 GeV [35] ) are used and the theoretical and experimental uncertainties are included; 
the knowledge of the precise value of this bound might be crucial to discriminate various 
scenarios of SUSY-breaking. We also discuss the implications of the LEP2 constraints 
on the MSSM Higgs sector and derive a lower bound on the value of tan/5, once various 
sources of uncertainty are taken into account. We then compare our results with those 
of previous analyses [3H1 EZl EH] which employ different computations of the Higgs boson 
masses, different codes for the determination of the (s)particle spectrum and which include 



only a subset of the errors discussed in the present analysis. 

The rest of the paper is organised as follows. In section 2, we describe the general 
structure of the spectrum generators, focusing on the aspects that are more critical to the 
mechanism of EWSB and the computation of the Higgs boson masses. In section 3 we com- 
pute the MSSM Higgs boson masses, compare the results from the three codes, and analyse 
the effect of various two-loop refinements, including the small effect of the 0(a.i,a. T ) and 
O(ol) corrections. We also discuss the theoretical uncertainties associated with the scheme 
dependence (comparing the DR and OS computations), the renormalisation scale depen- 
dence (variation with Mewsb), and the possible effect of the non-zero external momenta in 
the two-loop self-energies. In section 4 we study the impact of our results and the effects 
of the uncertainties on the upper bound on the lighter Higgs boson mass (and on its decay 
modes), as well as on the lower bound on the parameter tan/3, in various scenarios for the 
mechanism of SUSY breaking. We present our conclusions in section 5. Finally, in appendix 
A we provide for completeness the general formulae for the computation of the EWSB con- 
ditions and the Higgs masses, and in appendix B we provide explicit analytic formulae for 
the 0{a\,a T ) corrections to the MSSM Higgs boson masses and tadpoles. 

2 Determining the MSSM mass spectrum 

In this section, we describe the calculation of the MSSM mass spectrum implemented in 
the latest versions of three publicly available spectrum generators, i.e. SoftSusyl .8.7, 
SPheno2 .2.1 and SuSpect2 . 3. We start with the general algorithm used, then discuss 
some of the differences among the individual codes. The spectrum generator codes' manuals 
jHl EH 02] describe in detail the approximations used to perform their calculations. Here, 
we simply discuss the aspects of the calculation that are more critical to the Higgs boson 
masses. 

Assuming CP and R-parity conservation, the MSSM RGEs consist of some 110 coupled 
non-linear first-order homogeneous ordinary differential equations. The calculation of the 
MSSM mass spectrum consists of a two-boundary problem: we must solve these differential 
equations given boundary conditions at two values of the independent variable, the renor- 
malisation scale Q. Evolution of these RGEs is a standard problem, and there exist many 
techniques to perform the task [3D]- The high-scale (theoretical) boundary conditions are 
upon some of the parameters in the soft SUSY-breaking Lagrangian, and are specified by 
the model assumed for SUSY breaking (e.g., mSUGRA, AMSB, GMSB) and its parameters. 
They are often applied at the scale Mgut ~ O(10 16 ) GeV, at which the running electroweak 
gauge couplings unify: 

2i(M GUT )=3 2 (M GUT ), (1) 



where g\ = J5/3g' and g 2 = g with g and g' the couplings associated with the SU(2) L and 
U(l)y gauge groups of the MSSM, respectively. We do not enforce that the SU(3)c coupling 
constant g s = g^ unifies with g\ and g<i at the high scale; the few percent discrepancy 



from real unification is assumed to be accounted for by threshold corrections at Mgut- A 
case apart is the GMSB model, in which the high-scale boundary conditions are set at the 
messenger scale M mess . 

The weak-scale boundary conditions set the gauge and Yukawa couplings by matching 
the running MSSM parameters to the experimental data at some renormalisation scale, 
usually taken as Q = Mz- This step is quite involved, and requires subtracting the radiative 
corrections from the experimental data in order to arrive at the DR-renormalised MSSM 
parameters. Using the formulae of Ref. [I6 j and references therein, the MSSM DR gauge 
couplings g , g' , g% and the electroweak parameter v can be computed at Q — Mz from a set 
of four experimental input parameters. 

In the spirit of the SUSY Les Houches Accord (SLHA) for interfacing the spectrum 
generators with other computer codes [10], these input parameters can be chosen as: Gp, 
the Fermi constant determined from the muon decay; Mz, the pole mass of the Z boson; 
a era (M z ) MS , the five-flavour SM electromagnetic coupling at the scale Mz in the MS scheme; 
a s (Mz) MS , the five-flavour SM strong coupling at the scale M z in the MS scheme. In 
particular, the running parameters g , g' and v are connected to the running Z-boson mass 
mz by the relation: 

m z = Ml + ReIi T zz (Ml) = \^ + g")v\ (2) 

where U Z ' Z (M Z ) is the transverse part of the Z boson self-energy computed at a squared 
external momentum equal to the squared pole Z boson mass. 

The Yukawa couplings h u (u = u, c, t) for the up-type quarks, hd (d = d, s, b) for the 
down-type quarks and hi (£ = e, /i, r) for the leptons are determined from the corresponding 
running fermion masses as 

\[2m u V2m d sj^mi 

n u — : — a ' nd ~ a ' Ut - ~ a ' \ 6 ) 

v sinp v cosp v cosp 

where tan/3 = Vijv \ is the ratio of the vacuum expectation values (VEVs) of the two MSSM 
neutral Higgs fields, H± and H®, that also obey the relation v\ + v\ = v 2 . 

The running fermion masses m/ (with / = u, d, £) in eq. can be derived at the one-loop 
level from the corresponding pole masses Mf, through the relation 

m f = M f + Ef{Mf) (4) 

where S/(M/) is the one-loop fermion self-energy computed at an external momentum equal 
to the pole mass. In the case of the top quark, the self-energy includes also the leading two- 
loop standard QCD corrections 



m t 



M t + £ t (M t ) + (Am t ) 2 - lo °P<Q CD (5) 



where the precise expression for (Am t ) 2 1o °p>Q cd depends on the renormalisation scheme in 
which the parameters entering the one-loop self-energy E t are expressed. In the case of 



the bottom quark, the SLHA prescribes to take as input the SM running mass in the MS 
scheme, nib(nib) MS . Moreover, a "resummation" procedure is required (see e.g. Ref. |H]) in 
order to properly take into account the large QCD corrections, as well as the tan/3-enhanced 
SUSY corrections |%21 . to the relation between the input bottom mass and the corresponding 



MSSM, DR Yukawa coupling. We extract the latter, via eq. ©, from the MSSM, DR bottom 
mass nib, defined at the scale Q = Mz by the following matching condition: 

m b = m b (M z )™l SM = — — (6) 

where rf%b = m b (Mz)sM * s t ne SM, DR bottom mass, obtained by evolving mb{mb) MS up 
to the scale Q = M z with the appropriate RGE, in order to resum the QCD corrections, 
and then converting it to the DR scheme; A b = Ef,(mb) /nib accounts for the remaining 
non-gluonic corrections, some of which are enhanced by a factor tan/3. It has been shown 
pfTj that defining the running MSSM bottom mass as in eq. (JHJ) guarantees that the large 
threshold corrections of O (a s tan (3) n are included in nib to all orders in the perturbative 
expansion. In the case of the r lepton, the only tan/3-enhanced corrections to be included 
are those controlled by the electroweak gauge couplings, stemming from chargino-slepton 
loops. 

The computation of the MSSM mass spectrum is also complicated by the requirement 
that the electroweak symmetry is broken radiatively: at some weak scale M EWSB , the Higgs 
VEVs vi and v 2 can be computed by the minimisation of the MSSM effective potential, V e s ■ 
In principle, a set of high-energy boundary conditions is acceptable if, after RG evolution of 
the parameters down to M EW sb, it leads to the correct value of the squared running mass for 
the Z boson, as given in eq. (J2J). In practice, however, it is more convenient to assume that 
there is successful electroweak symmetry breaking, and trade two of the high-energy input 
parameters for v\ and t>2 (or, equivalently, for v and tan/3). The minimisation conditions 
of the effective potential then allow one to determine the weak-scale values of two MSSM 
parameters, usually chosen as the Higgs mass term in the superpotential, /x, and its soft 
SUSY-breaking counterpart, B, the latter being related to the squared running mass for the 
CP-odd Higgs boson A, niA, through the relation 

m 2 A = 2 B / sin2f3 . (7) 

The computation of /i and B beyond tree level requires the knowledge of the so-called 
tadpole corrections to the minimisation conditions of the effective potential. We provide 
the corresponding formulae in appendix A. At the one-loop level, Sof tSusy, SPheno and 
SuSpect use the complete calculation of the tadpoles given in Ref. [IB] . Among the two- 
loop corrections, the most relevant are those controlled by the top and bottom Yukawa 
couplings and by the strong gauge coupling. The three codes use the results of Ref. [2H] for 
the 0(a t a s + a>ba. s ) corrections, and those of Ref. [2B] for the 0{d^ + a t ab + ajj) corrections. 
Besides being required by consistency with the Higgs boson mass calculations, these two-loop 
corrections are relevant to stabilise the scale dependence of the MSSM parameters obtained 



from the minimisation conditions of the effective potential (see the discussion in Ref. 
Note that also for consistency reasons, the Higgs mass parameters m 2 H and rn 2 H2 which enter 
the minimisation conditions on V e s have to be evolved using the two-loop RGEs. 

After determining the full set of MSSM parameters, expressed in the DR renormalisation 
scheme at the scale M EW sb, the codes compute the spectrum of physical masses of the 
supersymmetric particles. The approximations employed in the computation of sfermion, 
gluino, chargino and neutralino physical masses differ from code to code, and we refer readers 
interested in the details to the respective manuals. 

On the other hand, due to extreme phenomenological relevance of the mass spectrum of 
the MSSM Higgs sector (as the failure of detecting a light Higgs boson already rules out a 
considerable fraction of the MSSM parameter space), it is mandatory to employ a precise 
calculation of the neutral Higgs boson masses that encompasses all the presently available 
radiative corrections. The general formulae for the computation of the Higgs masses are 
provided for completeness in appendix A. For the analysis presented in this paper, the 
spectrum generators SoftSusy, SPheno and SuSpect have been upgraded to include the 
state-of-the-art radiative corrections in the Higgs sector. 

In the newest versions Sof tSusyl .8.7, SPheno2 .2.1 and SuSpect2 . 3 the complete one- 
loop formulae, including the tadpole corrections, are taken from Ref. [ISj , which also provides 
formulae for the one-loop corrections to the charged Higgs boson mass. Concerning the two- 
loop part, the three codes employ the effective potential results of Ref. |22] for the 0(a t a s ) 
corrections, and those of Ref. |2H! for the 0(a 2 + a t ab + al) corrections. The formulae for the 
0(ab<y s ) corrections are obtained from those for the 0(at<y s ) corrections with appropriate 
replacements. The two-loop corrections are computed in the DR renormalisation scheme, i.e. 
they require that the one-loop part of the corrections is expressed in terms of running, DR 
renormalised MSSM parameters. This is a particularly convenient scheme in the situation 
in which the MSSM parameters are computed via RG evolution from a set of unified high- 
energy boundary conditions. 

The two-loop corrections to the Higgs boson masses controlled by the third-family 
Yukawa couplings can be completed by including the contribution of the r Yukawa cou- 
pling. The inclusion of the 0(a 2 ) corrections was discussed in Ref. 26J, and in appendix B 
we present a new computation of the remaining corrections of 0(ab<y T ). However, we found 
that the numerical impact of the r radiative corrections is very small in all the scenarios 
that we considered in our analysis. 

Other efforts to improve the two-loop calculation of the Higgs boson masses include an 
effective potential computation of the corrections to the lighter Higgs boson mass controlled 
by the electroweak couplings [3Qj, and a first attempt to go beyond the effective potential 
approximation by including the effect of non zero external momenta in the two-loop propaga- 
tors 43J. Two-loop corrections to the charged Higgs boson mass should also be calculated. 
We plan to include these results, as soon as they become available in a suitable form, in 
future updates of our spectrum generators. 

8 



It appears from the discussion above that the problem of determining the MSSM mass 
spectrum involves several mass scales, at which boundary conditions and EWSB conditions 
are imposed. As the computation of the radiative corrections at some scale requires knowl- 
edge about parameters that are determined at a different scale, an iterative procedure is 
adopted. The algorithm adopted by the codes can be summarised as follows: 

1. An initial guess for the MSSM parameters is taken at Q — Mz- 

2. The parameters are then evolved to Mqut (or M mess in the GMSB case), where the 
high scale boundary conditions are imposed. 

3. The new set of MSSM parameters is then evolved to M EW sb, where the EWSB con- 
straints are imposed. 

4. The MSSM pole mass spectrum is calculated. 

5. The parameters are then evolved to Q = Mz, where the Yukawa and gauge couplings 
are matched to empirical data. 

6. One then proceeds to step |2l continuing until the MSSM parameters converge upon 
stable values. 

Although the three codes employ the same general procedure in the determination of the 
MSSM mass spectrum, they differ in many details when it comes to the implementation of 
the single steps of the calculation. In particular, the differences that most affect the results 
for the Higgs boson masses are: 

- SPheno uses two-loop RGE for all the MSSM parameters, whereas SoftSusy and 
SuSpect use one-loop RGE for the soft SUSY-breaking masses and interaction terms 
of the sfermions (and two-loop RGE for the remaining MSSM parameters); 

- SoftSusy and SPheno include by default in the RGE the contributions of the Yukawa 
couplings for the three generations of fermions, whereas SuSpect works under the ap- 
proximation that the Yukawa couplings of the first two generations are zero. However, 
we found that the omission of the first two-family Yukawa couplings has always a very 
tiny impact on the Higgs boson masses, which is limited to the case of large tan /3; 

- SPheno defines the default EWSB scale M EW sb = sj m i± m t 2 i n terms of the physical 
stop masses, whereas SoftSusy and SuSpect define it in terms of the running stop 
masses (note, however, that all the codes also allow for an arbitrary choice of Mewsb); 

- when computing the running top mass at the scale Q = Mz according to eq. (jHJ), 
SPheno uses the running top mass itself in the one-loop part of the top self-energy, 
whereas SoftSusy and SuSpect use the pole top mass. This induces a difference in 
the formulae for the two-loop QCD corrections, which is properly taken into account 
by the codes; 

9 



- when computing the running bottom mass at the scale Q = Mz according to eq. (|5j). 
SPheno defines A& in terms of the full bottom self-energy, whereas SoftSusy and 
SuSpect include only the tan /3-enhanced contributions; 

- when computing the threshold corrections to the Yukawa couplings at the scale Q = 
Mz, SPheno uses the pole top mass in the Z boson self energy (required to obtain 
the running EWSB parameter v, see eq. (J2J)) whereas SoftSusy and SuSpect use the 
running top mass; 

- in the computation of the various threshold corrections at the scale Q = Mz, SoftSusy 
and SPheno use for all the sparticle masses and mixing angles the DR running values 
computed at Q = Mz, whereas SuSpect uses the values computed at Mewsb for the 
neutralino and chargino masses and mixing angles; 

- in the computation of the various radiative corrections at the scale Mewsb, SPheno 
and SuSpect compute v (Mewsb) through eq. (J2J), whereas SoftSusy takes the value 
computed at the scale Q = M z and evolves it up to M E wsb with the appropriate RGE. 

All of the differences listed above correspond to effects that are of higher order with 
respect to the accuracy required by the calculation of the Higgs boson masses. However, as 
will be discussed in the next section, they induce non negligible variations in the results for 
the Higgs boson masses, amounting to something less than a GeV for the lighter Higgs boson 
and a few GeV for the heavier Higgs bosons. For the time being, such variations must be 
taken as contributing to the uncertainties that still affect the calculation. Only when more 
precise computations of the radiative corrections become available will it be possible to fix 
the procedural ambiguities and reduce the uncertainty in the results. 

Before presenting our numerical results in the next section, we point out that the precise 
determination of the Higgs boson masses that we are discussing here does not only apply in 
constrained models such as mSUGRA, GMSB and AMSB. The neutral Higgs boson masses 
can also be determined in a general (unconstrained) MSSM, where the soft SUSY-breaking 
parameters are set by hand at the weak scale. The three spectrum generators have options 
in which such a procedure can be applied (although in the case of SoftSusy the procedure 
is not yet completely implemented). The programs can therefore also be viewed as DR 
Higgs boson mass calculators in a general MSSM framework, analogous to the FeynHiggs 
calculation in the On-Shell scheme 1 . At the end of the next section, we will present some 
illustrations of the calculation of the Higgs boson masses in the general or unconstrained 
MSSM. 



^dn fact, some of us are planning to include the part of the spectrum generators concerning the Higgs 
sector into the program HDECAY [23] which evaluates the Higgs boson decay widths and branching ratios in 
the SM and in the MSSM. 
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3 Precise results for the MSSM Higgs boson masses 

In this section we discuss in detail the results of the two-loop DR computation of the Higgs 
boson masses and of the EWSB conditions in the constrained MSSM (cMSSM), focusing on 
the mSUGRA, GMSB and AMSB models for SUSY breaking; we also display some results 
in the unconstrained MSSM. In a first part, we specify the physical scenarios in which we 
work. We then compare the results of the three codes, SoftSusy, SPheno and SuSpect, 
in the case of the cMSSM and discuss the impact of the various radiative corrections that 
affect the Higgs boson masses and the EWSB conditions. We subsequently discuss the 
various theoretical uncertainties which affect the Higgs boson mass calculations: the scheme 
dependence, the renormalisation scale dependence and the effect of the external momenta in 
the two-loop self-energies. We finally summarise the theoretical uncertainty on the Higgs 
boson mass determinations and discuss the impact of experimental errors on the SM input 
parameters on these calculations. 

3.1 The physical scenarios 

We will first work in the framework of constrained MSSM scenarios and illustrate our results 
in the case of the mSUGRA, GMSB and AMSB models for SUSY-breaking. For what 
concerns the choice of the high-energy boundary conditions, we restrict ourselves to six 
(out of ten) Snowmass Points and Slopes (SPS) [33]. Four of these points are mSUGRA 
ones, where the relevant input parameters are three universal soft SUSY-breaking terms, 
m o> m i/2 and Aq, the value of tan/5 (expressed in the DR scheme at the renormalisation 
scale Q = M z ) and the scale-invariant sign of //: 



SPS la 


: m 


= 100 GeV, 


mi/2 


= 250 GeV, 


A = -100 GeV, 


tan/? = 10 , 


/i>0, 


SPS 2 


m = 


= 1450 GeV, 


mi/2 


= 300 GeV, 


A = 0, 


tan/? = 10 , 


/i>0, 


SPS 4 


m 


= 400 GeV, 


mi/2 


= 300 GeV, 


A = 0, 


tan/? = 50 , 


/i>0, 


SPS 5 


m 


= 150 GeV, 


mi/2 


= 300 GeV, 


A) = -lTeV, 


tan /3 — 5 , 


n> 0. 



The first choice, SPSla, corresponds to a "standard" mSUGRA scenario; the second 
choice, SPS2, is called a "focus point" scenario and has interesting implications for dark 
matter abundance; SPS4 is characterised by a large value of tan j3 while SPS5 is characterised 
by a large value of the stop mixing parameter and leads to a light stop squark. The focus point 
and the large tan f3 scenario are known to be the most sensitive to the various approximations 
made by the spectrum generators [11J. 

We will also choose one point in the GMSB scenario, where the input parameters are the 
SUSY-breaking scale A, the messenger scale M mess , the messenger index N mesB , tan/? and 
the sign of fi: the SPS8 point in which the lightest SUSY particle is the lightest neutralino 

SPS 8 : A = 100 TeV, M mcss = 200 TeV, N mcss = 1 , tan/? = 15 , // > . 
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The last point, denoted SPS9, is for an AMSB scenario described by a large gravitino 
mass 777.3/2, a common mass term for the scalars m , tan/? and the sign of \x: 

SPS 9 : m 3/2 = 60 TeV, m Q = 450 GeV, tan/3 = 10 , pt > . 

We remark in passing that in Ref. [31] the SPS scenarios are defined in terms of the 
low-energy MSSM DR parameters, as computed by the code IsaJet7.58 at the weak scale 
Mewsb- Thus, calling "SPS scenarios" the above choices of high-energy boundary conditions 
always involves a slight abuse of language. 

In a next step, we will work in the case of the unconstrained MSSM. For the latter 
scenario, we take the so-called "phenomenological MSSM" (pMSSM) [2j, in which some 
constraints have been imposed (such as CP conservation, flavour diagonal sfermion mass 
and coupling matrices and universality of the first and second generations). These con- 
straints ensure that an appreciable fraction of parameter space of the pMSSM has viable 
phenomenology. The model involves 22 free parameters in addition to those of the SM: 

- the parameter tan/3; 

- the two soft SUSY-breaking Higgs mass parameters m 2 Hi and m 2 H2 , which can be traded 
against Ma and fi by enforcing the EWSB conditions; 

- the three gaugino mass parameters M% , M 2 and M 3 ; 

- the diagonal sfermion mass parameters m? : five for the third generation sfermions 
and five others for the first/second generation sfermions; 

- the trilinear sfermion couplings Af. three for the third generation sfermions and three 
others for the first/second generation sfermions. 

Fortunately, most of these parameters have only a marginal impact on the Higgs boson 
masses, a fact which considerably simplifies the analysis. A very important parameter in 
the context of the radiative corrections to the MSSM Higgs boson masses is the mixing 
parameter in the stop sector, X t = A t — \x cot (3. For our numerical illustrations, we will 
consider three pMSSM scenarios, in which the SUSY spectrum is rather heavy: the soft 
SUSY-breaking masses of the sfermions (all set equal to a common mass Ms), the SU(3) 
gaugino mass M 3 (Mi and M 2 being related to the latter through the usual GUT relation), 
the pseudoscalar Higgs boson mass Ma and the higgsino mass parameter /i are set to 1 TeV; 
for tan/3 we choose a moderate value, tan/3 = 10; furthermore, we consider three cases for 
the mixing X t in the stop sector: no mixing, X t = 0, typical mixing, X t = Ms, and large 
mixing, X t = \/QMs (the other sfermion trilinear couplings will be also set to 1 TeV). All 
of the MSSM input parameters are taken as DR running quantities computed at the scale 
Q = 1 TeV, apart from tan/3 which is computed at Q = Mz, and Ma which denotes the 
physical mass. Thus we will have three pMSSM points: 

H = 1 TeV, tan/3 = 10, X t = , 
(j. = 1 TeV, tan/3 = 10, X t = M s , 
fj, = 1 TeV, tan/3 = 10, X t = V6M S . 
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pMSSMl : 


M s = A b:T 


= M 3 


= M A 


pMSSM2 : 


M s = A b , T 


= M 3 


= M A 


pMSSM3 : 


M s = A b , T 


= M 3 


= M A 



The first and third of these scenarios are rather close to those which have been used 
as benchmark points in the interpretation in the unconstrained MSSM of the Higgs boson 
searches by the LEP2 collaborations [35] . 

Finally, the SM input parameters at the weak scale are fixed according to the SLHA 
prescription. In particular, we take for the electroweak and strong parameters [46J: 

G F = 1.16639 1(T 5 GeV" 2 , M z = 91.1876 GeV, 

a~^{M z ) m = 127.934 ± 0.027, a s (M z ) m = 0.1172 ± 0.002, (8) 

and for the third-generation fermion masses the values [331 I46j : 

M t = 178.0 ± 4.3 GeV, m b (m b ) m = 4.25 ± 0.25 GeV, M T = 1.777 GeV. (9) 

In most of our discussion we take the SM input parameters to be equal to their central 
values (the errors on Gf, Mz and M T are indeed very small and we omit them in the 
equations above). However, in section I3~B1 we discuss the uncertainty in the determination 
of the Higgs masses arising from the experimental errors on some of these parameters. 

3.2 Determination of the neutral Higgs masses in the cMSSM 

3.2.1 Comparing the results of the three codes 

We start our discussion by comparing the values of the neutral Higgs boson masses and of 
the superpotential parameter /i in the six SPS scenarios discussed above, as they result from 
the computations of SoftSusy, SPheno and SuSpect at the default renormalisation scale 
Mewsb = v/wij Tf% , once all of the two-loop radiative corrections are implemented. The 
corresponding values for the physical masses of the CP-even Higgs bosons, Mh and Mh, 
for the physical mass of the CP-odd Higgs boson, M^, and for the parameter fj, (the latter 
interpreted as a DR running parameter computed at Mewsb) are given in the tables HJ El 
and|H respectively. 

Concerning the lighter CP-even Higgs boson mass, M^, it can be seen that the agreement 
between the three codes is very good, the discrepancies being contained in a half GeV. In the 
cases of the heavier Higgs boson masses and of /i, which are more sensitive than Mh to small 
variations in the RG evolution of the soft SUSY-breaking parameters, the discrepancies 
amount to at most few GeV, generally below the 1% level. We find this agreement very 
satisfactory, taking into account the fact that the results are obtained with three independent 
codes that - although based on the same set of formulae for the corrections to the Higgs 
boson masses and the EWSB conditions - differ in many details of the calculation, as outlined 
in the previous section. 

We have checked that, if we force the three codes to use the same computation of the 
threshold corrections to the gauge and Yukawa couplings and the same set of RGE, the 
residual discrepancies in the results for the neutral Higgs boson masses and \i become neg- 
ligible. However, we stress again that the differences between the three codes are a matter 
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of choice, because they all correspond to effects that are of higher order with respect to the 
accuracy required by the calculation. The spread in the values of the Higgs boson masses 
and \x resulting from tables EE1 has to be considered as part of the theoretical uncertainty, 
and will be reduced only when more refined calculations of the radiative corrections become 
available. 



Code 


SPSla 


SPS2 


SPS4 


SPS5 


SPS8 


SPS9 


SoftSusy 

SPheno 

SuSpect 


112.1 
112.2 
112.1 


116.8 
117.1 
116.8 


114.1 
114.3 
114.1 


116.3 
116.5 
116.1 


115.4 
115.8 
115.5 


117.4 
117.8 
117.5 



Table 1: The lighter CP-even Higgs boson mass, Mh, in the six SPS scenarios, as computed 
by Sof tSusyl .8.7, SPheno2 .2.1 and SuSpect2 . 3. The SM input parameters are chosen as 
in eqs. ©-(0). 



Code 


SPSla 


SPS2 


SPS4 


SPS5 


SPS8 


SPS9 


SoftSusy 

SPheno 

SuSpect 


406.5 
406.0 
406.5 


1553.0 

1554.6 
1552.1 


355.8 
360.5 
355.3 


686.8 
686.5 
686.9 


550.4 

552.4 
550.6 


1056.9 
1051.1 
1056.6 



Table 2: Same as tabled for the heavier CP-even Higgs boson mass, Mh- 



Code 


SPSla 


SPS2 


SPS4 


SPS5 


SPS8 


SPS9 


SoftSusy 

SPheno 

SuSpect 


406.2 
405.7 
406.1 


1552.9 

1554.5 
1552.0 


355.8 
360.5 
355.3 


687.0 
686.9 
687.2 


550.1 

552.1 
550.3 


1056.8 
1051.0 
1056.5 



Table 3: Same as table [T] for the CP-odd Higgs boson mass, Ma- 



Code 


SPSla 


SPS2 


SPS4 


SPS5 


SPS8 


SPS9 


SoftSusy 

SPheno 

SuSpect 


364.8 
364.3 
364.7 


586.5 
588.2 
583.6 


413.8 

414.7 
413.6 


631.2 
631.2 
631.3 


440.1 
442.2 
440.3 


1011.8 
1005.9 
1011.1 



Table 4: Same as table Q for the superpotential Higgs mass parameter /i, expressed in the 
DR renormalisation scheme at the scale Mewsb- 
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3.2.2 Impact of the radiative corrections 

We now discuss the importance of the various radiative corrections to the Higgs boson masses 
and tadpoles that are taken into account by the codes SoftSusy, SPheno and SuSpect. 
Tables El and El show the variations of Mh and Mh, respectively, in the six SPS scenarios, as 
a result of different approximations for the radiative corrections. The first three rows of each 
table contain the values of the masses as obtained by SuSpect employing, respectively, the 
tree-level, one-loop and two-loop formulae for the Higgs masses and EWSB conditions. The 
variations shown in the second part of the tables result from progressively switching on the 
various two-loop corrections to the Higgs boson masses and EWSB conditions. Explicitly, 
we examine the two loop 0(a t a s + a b a s ) corrections (first line), the 0(a% + a t a b + of) 
corrections (second line) and the newly calculated 0(a^ + a b a T ) corrections (third line). 



Approximation 


SPSla 


SPS2 


SPS4 


SPS5 


SPS8 


SPS9 


Tree-level 


88.6 


87.8 


90.2 


82.9 


89.2 


88.2 


One-loop 


109.7 


112.3 


111.1 


113.0 


111.3 


112.6 


Two-loop 


112.1 


116.8 


114.1 


116.1 


115.5 


117.5 


0(a t a s + a b a s ) 


+3.1 


+5.7 


+3.8 


+3.0 


+5.3 


+5.8 


0[a\ + a t a b + af) 


-0.6 


-1.1 


-0.8 


+0.2 


-1.1 


-1.0 


0{a 2 T + aba T ) 


< 10~ 3 


< 10~ 3 


< 10- 3 


< io- 3 


< 10" 3 


< 10~ 3 



Table 5: The lighter CP-even Higgs boson mass, Mh, in the six SPS scenarios, as computed 
by SuSpect under different approximations for the radiative corrections. The first three rows 
contain the mass (in GeV) computed at tree-level, one- and two-loop, respectively; the last 
three rows contain the shifts (in GeV) due to the different two-loop contributions. 



Approximation 


SPSla 


SPS2 


SPS4 


SPS5 


SPS8 


SPS9 


Tree-level 


402.5 


1542.7 


356.3 


689.6 


525.7 


1043.1 


One-loop 


406.8 


1551.3 


356.0 


688.2 


549.0 


1056.1 


Two-loop 


406.5 


1552.1 


355.3 


686.9 


550.6 


1056.6 


0(a t a s + a b a s ) 


-0.1 


+ 1.1 


-0.1 


-0.9 


+2.5 


+0.9 


0(af + a t a b + 0%) 


-0.3 


-0.3 


-0.5 


-0.3 


-0.9 


-0.4 


0{o? T + a b a T ) 


< 10~ 3 


< 10~ 3 


0.01 


< 10~ 3 


< 10~ 3 


< 10~ 3 



Table 6: Same as table El for the heavier CP-even Higgs boson mass, Mh 

Comparing the first three rows of table El it can be seen that the one-loop corrections 
to Mh are positive and of the order of 20 GeV, while the two-loop corrections to Mh are 
also positive and of the order of 2-5 GeV. This is to be contrasted with the case of the 
OS calculations (see e.g. Ref. [23) 5 where both the one-loop and two-loop corrections are 
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usually larger, and the two-loop corrections are negative. However, as will be discussed 
in section 13.3.11 the two-loop results of the OS and DR calculations agree within 2 GeV 
in general. From the first three rows of table El instead, it can be seen that the one-loop 
corrections to Mh amount in general to few GeV, and their impact is usually of the order of 
1% (apart from the scenario SPS8, where they reach 4%). The two-loop corrections are also 
small, being of the order of one GeV or less. This is due to the fact that, in the considered 
scenarios, the heavier Higgs boson masses are large already at the tree-level, so that the 
inclusion of radiative corrections can have only a small effect. 

Further information on the relative size of the various two-loop corrections to the Higgs 
boson masses can be obtained from the lower parts of tables El and El In the case of Mh, 
we see that the 0(a t a s + ab<y s ) corrections, controlled by the strong gauge coupling, are 
dominant (amounting to 3-6 GeV) and positive, whereas the 0{a"l + a t ab + a%) corrections, 
involving only the top and bottom Yukawa coupling, are smaller (of the order of one GeV) 
and usually negative (apart from the scenario SPS5). This is again to be contrasted with the 
case of the OS calculations, where the 0(ottot s ) corrections, which are always leading at two- 
loop order, are negative, and the subleading 0(af) corrections are positive (in both OS and 
DR calculations, the corrections controlled by the bottom Yukawa coupling are negligible 
for small or moderate values of tan/3). In the case of Mh, instead, the trend is less definite: 
the 0(a t a s + af,a s ) corrections can be either positive or negative in the different scenarios, 
and they are not always larger than the 0(a^ + a t ab + af) ones. 

Concerning the corrections controlled by the r Yukawa coupling, it can be seen that their 
effect is usually negligible, amounting to less than a MeV for both Mh and Mh- Only in 
the scenario SPS4, in which the r Yukawa coupling is enhanced by a large value of tan/?, 
a very tiny (10 MeV) effect is visible in the heavier Higgs boson mass. This suppression of 
the r corrections even in the case of large tan /3 is explained by the smallness of m T (even in 
comparison with m&) , by the absence of colour enhancements, and by the fact that the only 
tan /3-enhanced threshold corrections to the relation between h T and m T are those controlled 
by the electroweak gauge couplings. 

One has to bear in mind, however, that all of the above considerations on the relevance of 
the various two-loop corrections to Mh and Mh depend to some extent on the choice of the 
electroweak symmetry breaking scale, taken by default as Mewsb = ^/m^ m 4 ~ 2 . For example, 
as will be clear from the figures in section 13.3.21 with the choice M EW sb = M z the two-loop 
corrections to Mh become smaller (1-2 GeV) and negative, while the two-loop corrections to 
Mh amount to several GeV and they are positive. Also, in the scenario SPS4 the corrections 
to Mh controlled by the r Yukawa coupling get somewhat larger (~ 0.2 GeV). It is thus 
clear that - although it is always possible to choose Mewsb in such a way that some of 
the two-loop corrections are small - only the proper inclusion of all of the leading two-loop 
corrections to the Higgs boson masses allows us to obtain reliable results that do not depend 
strongly on a preconceived choice of the renormalisation scale. 
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3.2.3 Impact of some higher— order effects 

To understand the importance of several higher-order effects that are taken into account by 
SoftSusy, SPheno and SuSpect, it may be useful to discuss their individual impact on the 
numerical values of the Higgs masses. 



Approximation 


SPSla 


SPS2 


SPS4 


SPS5 


SPS8 


SPS9 


Default 


112.17 


117.06 


114.34 


116.46 


115.80 


117.82 


2-loop sfermion RGE 


0.07 


0.15 


0.04 


0.03 


0.03 


0.03 


No Am t 21oop 


0.26 


0.33 


0.28 


0.35 


0.35 


0.29 


No Amf w 


0.14 


0.40 


0.50 


0.15 


0.44 


0.29 


No m& resummation 


< 10~ 3 


< 10" 3 


0.06 


< 10" 3 


< 10~ 3 


< 10~ 3 


Pole Higgs masses in 1-loop 


0.02 


0.02 


0.01 


0.01 


0.02 


0.02 



Table 7: Lighter CP-even Higgs mass, Mh, in the six SPS scenarios, as computed by SPheno 
under different approximations for the higher-order effects. The first row contains the default 
values in GeV, the other rows contain the shifts (in GeV) due to the approximations. 



Approximation 


SPSla 


SPS2 


SPS4 


SPS5 


SPS8 


SPS9 


Default 


406.0 


1554.6 


360.5 


686.5 


552.4 


1051.1 


2-loop sfermion RGE 


1.1 


6.8 


0.3 


-1.3 


0.7 


1.5 


No Amf oop 


1.1 


6.7 


4.5 


-1.4 


2.6 


-2.1 


No Amf w 


0.6 


8.2 


8.0 


-0.6 


3.2 


-2.2 


No rrib resummation 


0.2 


0.3 


31.0 


0.1 


0.3 


<io- 2 


Pole Higgs masses in 1-loop 


0.1 


0.1 


0.5 


0.1 


0.3 


0.1 



Table 8: Same as table [H for the heavier CP-even Higgs boson mass, Mh- 



Tables and |H1 show the variations of Mh and Mh, respectively, in the six SPS scenarios, 
as a result of different approximations for the higher-order corrections that affect the RGE 
evolution of the MSSM parameters. The first row of each table contains the default values 
of the masses, as obtained with the public version of SPheno2.2.1. In the first row of 
the second part of the table, we show the effect (in GeV) from switching off the two-loop 
terms in the RGEs for the sfermion soft SUSY-breaking parameters; in the second row, from 
switching off the two-loop QCD corrections in the computation of the top Yukawa coupling 
h t ; in the third, from switching off the one-loop electroweak corrections in the computation 
of ht] in the fourth, from defining the running bottom mass as nib — ^b (1 + A^) instead of 
adopting the "resummed" formula of eq. (jfi|); in the fifth, from using the physical values of 
the Higgs masses, instead of the running values, as input parameters for the computation of 
the one-loop corrections. 
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We see that the effect of these approximations on Mh is in general moderate, amounting 
at most to some tenths of a GeV, whereas their effect on Mjj can amount to many GeV. 
This is due to the fact that for Ma 3> Mz, as is always the case in the SPS scenarios, the 
tree-level value of Mh is essentially proportional to mz, and the dependence on the MSSM 
parameters enters only through the radiative corrections. On the other hand, the heavier 
Higgs boson masses depend already at tree-level on the soft SUSY-breaking parameters m 2 Hl 
and m 2 H2 , and the RG-evolution of the latter is very sensitive to the precise value of all of 
the other parameters, especially h t . We also see that, as expected, the precise definition of 
the bottom Yukawa coupling has a sizeable effect on the Higgs boson masses only in the 
scenarios with large tan j3. Finally, we see that the choice of the renormalisation scheme for 
the values of the Higgs boson masses entering the one-loop corrections also induces a small 
variation in the final results for the physical Higgs boson masses. We note, however, that 
this variation amounts to a two-loop effect controlled by the electroweak gauge couplings, 
thus it is of higher order with respect of the accuracy of our calculation. 

3.3 Estimate of the theoretical uncertainties in the cMSSM 

In this subsection, we will discuss the shifts on the MSSM neutral Higgs boson masses 
induced by the choice of a different renormalisation scheme and by the variation of the 
renormalisation scale. We will also try to estimate the impact of the approximation of zero 
external momentum in the two-loop corrections. These shifts indicate the size of the higher 
order radiative corrections that are left uncomputed, and thus give a rough estimate of the 
theoretical uncertainties in the calculation of the masses. 



3.3.1 Comparing the DR and OS calculations 



The first interesting comparison to make is between the results of our DR calculation of the 
CP-even Higgs boson masses and those of the OS calculation implemented in FeynHiggs. 
The latter is based on the one-loop formulae of Refs. (T3J E2] and the two-loop formulae 
of Ref. |Uj (for the 0(a t a s ) part) and Refs. J2H 123 H3 (for the 0(a b a s ), 0(a 2 ) and 
0(a t ab+al) parts, respectively). In the computation of FeynHiggs, the one-loop corrections 
are expressed in terms of physical masses and mixing angles, and appropriate counterterm 
contributions are inserted in the two-loop corrections for consistency. 

As FeynHiggs does not perform the evolution of the MSSM parameters from the high- 
energy input scale to Mewsb, in order to obtain the Higgs boson masses in the SPS scenarios 
we have to provide it with the input parameters at the weak scale, as they are computed 
by one of the RGE codes. In particular, the top quark mass must be the physical one, 
M t ; the parameters \i and tan/3 must be the running ones, computed at the scale Mewsb; 
the squark masses and mixings can be either physical or running (in the latter case they 
are converted by FeynHiggs into the physical ones by applying the appropriate one-loop 
corrections); the definition of the gaugino mass parameters is irrelevant at the perturbative 
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order we are interested in; finally, the A boson mass, which enters the tree-level mass matrix 
of the CP-even Higgs bosons, must be the physical one, computed at the two-loop accuracy. 

Table El shows the values of M^ and Mh as computed by FeynHiggs (version 1.5.1) in 
the six SPS scenarios, with the weak scale input parameters provided by SuSpect. The 
corresponding results from the pure SuSpect calculation are also shown for comparison. It 
can be seen from tableElthat the values of Mh resulting from the OS calculation of FeynHiggs 
are consistently a couple of GeV larger than the ones provided by the DR calculation. As the 
DR and OS computations are equivalent up to the two-loop 0(a t a s + an,a s + a>1 + a. t otb + a l) 
accuracy, the discrepancies must be due to terms that are formally of higher order, i.e. 
two-loop terms controlled by the electroweak gauge couplings and three-loop terms, among 
which the most important are controlled by the top Yukawa and strong gauge couplings. 
Concerning the heavier Higgs boson mass, Mh, it can be seen in table |H1 that the results of 
FeynHiggs are in excellent agreement with those of SuSpect. This fact is not surprising, 
because for Ma ^> M z , as is always the case in the SPS scenarios, M H and Ma are very 
strictly correlated, and in our comparison FeynHiggs takes as input the value of Ma as 
computed by SuSpect. 



Mass 


Code 


SPSla 


SPS2 


SPS4 


SPS5 


SPS8 


SPS9 


M h 


SuSpect 
FeynHiggs 


112.1 
113.8 


116.8 
118.3 


114.1 
116.1 


116.1 
118.5 


115.5 
117.3 


117.5 
118.3 


M H 


SuSpect 
FeynHiggs 


406.5 
406.5 


1552.1 
1552.0 


355.3 

354.8 


686.9 
686.5 


550.6 
550.6 


1056.6 
1056.7 



Table 9: The lighter and heavier CP-even Higgs boson masses, Mh and Mh, in the six 
SPS scenarios, as computed by SuSpect and FeynHiggs. The weak scale input parameters, 
including the physical mass Ma, are taken from the output of SuSpect. 



3.3.2 Renormalisation scale dependence 

Another measure of the effect of the higher orders consists in studying the numerical depen- 
dence of the results for the physical Higgs masses on the renormalisation scale M EWSB at 
which the effective potential is minimised and the radiatively corrected masses are computed. 
In the ideal case of an all-orders calculation, the physical observables should not depend on 
the choice of the scale. The residual scale dependence still present in the real case can be 
taken as a rough estimate of the magnitude of the corrections that are left uncomputed. 
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Figure 1: The lighter CP-even Higgs boson mass, M^, as computed by SoftSusy in the 
six SPS scenarios, as a function of the minimisation scale of the effective potential, Mewsb- 
The dotted lines correspond to the one-loop computation and the solid lines to the two-loop 
computation. 
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Figure 2: Same as figure Q for the heavier CP-even Higgs boson mass, Mjj- 



21 



Figs. ^ and El show the renormalisation scale dependence of the CP-even Higgs boson 
masses M h and M H , respectively, as computed by SoftSusy in the six SPS scenarios dis- 
cussed in section 13.11 The scale dependence of Ma is very similar to that of Mh and is thus 
omitted for brevity. We have also checked that the same features can be reproduced with 
SPheno and SuSpect. The dotted line in each plot corresponds to the one-loop computation 
of the relevant mass, whereas the solid line corresponds to the two-loop computation. The 
scale Mewsb is varied from a minimum value of Mz to a maximum corresponding to twice 
the stop mass scale Jrni x m^ 2 . 

It can be seen from fig. ^ that the one-loop results for the lighter Higgs boson mass 
Mh show a sizeable scale dependence, varying by nearly 10 GeV in the considered range 
of M EWSB . This dependence is essentially driven by the variation of the running top mass, 
which decreases as the renormalisation scale increases (we recall that the leading one-loop 
corrections to M 2 are proportional to m\jv 2 \ On the other hand, it appears that the 
inclusion of the two-loop corrections significantly improves the scale dependence of Mh, 
compensating for the dependence on m t and leaving a residual variation of 2-3 GeV in the 
considered range of M EW sb- 

Fig. El shows that a similar situation occurs in the case of the heavier Higgs boson mass 
Mfj, with the inclusion of the two-loop corrections clearly improving the renormalisation 
scale dependence of the results. One exception is given by the plot corresponding to the 
scenario SPS4, where the curves show a spike and the two-loop corrections do not improve 
the scale dependence. This anomalous behaviour is due to the fact that the codes use the 
running values of the Higgs boson masses in the computation of the one-loop corrections. 
For large values of tan/3 (e.g. tan/3 = 50 in SPS4) the running masses for the heavier Higgs 
bosons vary by hundreds of GeV as the renormalisation scale varies in the considered range. 
Thus, for some values of the renormalisation scale, the loop integrals involving the heavier 
Higgs boson masses encounter real-particle thresholds which give rise to the characteristic 
spike behaviour. We have checked that, if we use in the one-loop corrections the physical 
values of the Higgs boson masses, which vary only of a few GeV in the considered range 
of Mewsb, the spike disappears and the plot for the SPS4 scenario looks qualitatively very 
similar to the ones for the other scenarios (the same happens for another large-tan /3 scenario, 
SPS lb, which we did not include in our analysis). We recall, however, that the choice of 
whether to use the running or physical Higgs masses in the one-loop corrections amounts to 
a higher-order effect with respect to the accuracy of our calculation. 

3.3.3 The external momentum dependence 

To simplify the calculation of the two-loop radiative corrections to the Higgs boson masses, 
the squared external momentum p 2 in the two-loop part of the self energies has been set to 
zero, rather than to the squared mass of the corresponding Higgs boson. This approximation 
can be considered as a source of theoretical uncertainty. The actual size of the shift on the 
Higgs boson masses induced by this approximation can only be obtained when the exact two- 
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loop computation, retaining the external momentum dependence, is performed (see Ref. 
for a first attempt). However, it can be roughly estimated by considering the error that the 
zero-momentum approximation induces in the one-loop corrections, and assuming that the 
error induced in the two-loop corrections will have the same relative size (a similar analysis 
was performed in Ref. |33J). 



Mass 


Approximation 


SPSla 


SPS2 


SPS4 


SPS5 


SPS8 


SPS9 


M h 


p 2 = in 1-loop 
Estimated 2-loop effect 


+0.9 
0.1 


+2.0 
0.4 


+ 1.3 
0.2 


+1.0 
0.1 


+2.0 
0.3 


+1.7 
0.3 


M H 


p 2 = in 1-loop 
Estimated 2-loop effect 


-0.3 
0.02 


+0.8 
0.1 


+4.4 
0.7 


+0.02 
0.02 


-0.1 
<0.01 


-0.9 
0.04 



Table 10: Effect of the zero-momentum approximation in the one-loop corrections to the 
CP-even Higgs boson masses, Mh and Mh, in the six SPS scenarios, and an estimate of the 
corresponding two-loop effect. All the numbers are in GeV. 

In table EH we display the effect of the approximation of zero external momentum on 
the lighter and heavier CP-even Higgs boson masses, as computed by SuSpect in the six 
SPS scenarios. For each mass, the first row contains the shift (in GeV) resulting from the 
approximation p 2 = in the one-loop part of the corrections, and the second row contains the 
estimate for the size of the same effect in the two-loop case. The latter is obtained by scaling 
the absolute value of the shift induced in the one-loop corrections by the ratio between the 
size of the two-loop zero-momentum corrections and that of the one-loop zero-momentum 
corrections in the considered scenario. As can be seen, the effect of the zero-momentum 
approximation on the one-loop calculation of Mh is of the order of 1-2 GeV. When rescaled 
to account for the relative size of the two-loop to the one-loop corrections, this results in a 
shift of less than half GeV in all of the considered scenarios. 

On the other hand, we see from table E3 that in all scenarios but SPS4 the effect of the 
zero-momentum approximation on the one-loop computation of Mh is of the order of a few 
tenths of GeV. As a consequence, one can predict that in those scenarios the effect of the 
same approximation on the two-loop corrections is negligible. This is due to the fact that, 
as mentioned in section I3.2.2| the heavier Higgs boson mass is dominated by the tree-level 
value, and the impact of the one-loop radiative corrections is anyway small. In the case 
of the large-tan /3 scenario SPS4, however, there are non-negligible momentum-dependent 
contributions from the one-loop diagrams with bottom quarks and tau leptons, inducing a 
shift of more than 4 GeV in Mh, and the effect of the zero-momentum approximation on 
the two-loop corrections can be predicted to get close to one GeV. 
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3.4 The case of the pMSSM 

We now discuss the case of the phenomenological MSSM defined in section ETT1 In principle, 
since there are 22 free parameters in the model, the phenomenological analyses should be 
rather complicated to carry out. However, as mentioned previously, only a small subset 
of parameters plays a significant role in the evaluation of the Higgs boson masses. Before 
presenting our numerical results, we briefly summarise the role and the main effects of the 
various MSSM parameters. 

At the tree level, the Higgs sector of the pMSSM can be described by two input parameters 
in addition to the SM ones; these parameters are in general taken to be the mass of the 
pseudoscalar Higgs boson Ma and tan/3. In this parameterisation, the tree-level CP-even 
Higgs boson masses are given by 



M k/H - -L 



M\ + M 2 Z T \j(M\ + M z ) 2 - 4 cos 2 2pM\M 2 z 



1/2 



(10) 



The charged Higgs boson mass is given at tree-level by M#± = \jM\ + Myy, while the 
angle a in the CP-even Higgs sector is given by tan 2a = tan 2/3 {M\ + M Z )/{M\ — Mf), 
with the constraint — | < a < 0. 

The maximal value of the lighter CP-even Higgs boson mass at the tree level, Mh = 
cos 2/3 M z , is obtained in the decoupling regime, Ma 3> M z (note that, in this case, the 
lighter Higgs boson of the MSSM will have exactly the same properties as the SM Higgs 
boson), and is saturated for large values of tan/3. In this limit, the mass of the heavier 
CP-even Higgs boson is simply Mh = Ma- 

The leading one-loop radiative corrections to M 2 are controlled by the top and bottom 
Yukawa couplings, and in the decoupling regime read 

3m\ / Ml , Xf Xf \ 3mf X* 



In -11- + JUL —t o o (-\-\\ 

' 2n 2 v 2 \ m 2 2M| 12MIJ 2tt 2 v 2 12MI ' l ; 

where Ms is a common soft SUSY-breaking mass term for the third-generation squarks and 
X t fi are given in terms of the stop/sbottom trilinear couplings, \x and tan/3, by: 

X t = A t — /i cot /3 and X b = A b — /itan/3 . (12) 

For M A 3> M z and tan/3 3> 1, one obtains M\ = M z + e for the lighter CP-even 
Higgs boson, while for the heavier Higgs boson one still has Mh = Ma- The choice of 
the renormalisation scheme for the parameters entering eq. ()11|) becomes relevant when the 
two-loop corrections are included; we recall that in our calculation we express the one-loop 
corrections in terms of DR parameters. 

Because of the quartic dependence on m t , the corrections controlled by the top Yukawa 
coupling are the leading ones. These corrections are larger when the logarithm in the first 
term of eq. (|TT|) is larger, i.e. for large M s values (corresponding to large stop masses). In 
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addition, the top quark corrections are maximal in the so-called M™ ax scenario |45j . where 
the trilinear stop coupling in the DR scheme is such that X t ~ v6Ms. 

The corrections controlled by the bottom Yukawa coupling are in general strongly sup- 
pressed with respect to those controlled by the top Yukawa coupling, due to the overall 
factor m\. However, in the last term of eq. (|11|). proportional to X£, this suppression can 
be compensated by a large value of the product /itan/3, providing a non-negligible negative 
correction to M%. The choice of the values for the remaining soft SUSY-breaking parameters 
does not have a very large impact on the one-loop corrections, and in the DR calculation 
the two-loop corrections, although numerically significant in the determination of the precise 
value of the lighter Higgs boson mass, do not substantially alter the picture. 
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Figure 3: The lighter MSSM Higgs boson mass as a function of X t in the DR scheme for 
tan/3 = 10 and M s = Ma = 1 TeV with M t = 178 GeV. The full and dashed lines correspond, 
respectively, to the two-loop and one-loop corrected masses as calculated with SuSpect, 
while the dotted line corresponds to the two-loop Mh value obtained with FeynHiggs. 

The above features are exemplified in fig. EJ where the lighter Higgs boson mass is dis- 
played as a function of the DR parameter X t , for M t = 178. In the figure, the MSSM 
parameters are set to those of the three pMSSM points introduced in section 13. 1| in par- 
ticular, the physical pseudoscalar mass Ma and the third-generation soft SUSY-breaking 
scalar masses Ms (the latter computed at the renormalisation scale Q = 1 TeV) are set to 
1 TeV, while tan/3 is fixed to tan/3 = 10 at Q = Mz- The dashed curve for the one-loop 
corrections, and the full curve for the two-loop corrections in the DR scheme, have been 
obtained using the program SuSpect. As one can see, the lighter Higgs boson mass Mh 
has a local minimum for zero stop mixing, and it increases with \X t \ until it reaches a local 
maximum at the point \X t \ = yQMs ~ 2.45 TeV, where it starts to decrease again. Note 
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that the two-loop corrections amount to several GeV and, in this scenario, decrease in size 
when X t takes on large and positive values. 

The dotted curve in fig. |3] is obtained with the program FeynHiggs, using the MSSM 
input parameters as they come from SuSpect. In the OS renormalisation scheme adopted by 
FeynHiggs, the maximal Mh value is obtained for X® ~ 2M® S , where X® s and M® s are 
the unphysical parameters obtained by rotating the diagonal matrix of the OS stop masses 
by the OS mixing angle (see e.g. Ref. [32j f° r a discussion). On the other hand, as we are 
plotting Mh as a function of the DR parameter X t for both the SuSpect and FeynHiggs 
computations, we find the maximum value of Mh roughly at the same place. Comparing the 
solid and dotted lines, it can be seen that for small to moderate values of X t the results of the 
OS calculation of FeynHiggs are a couple of GeV higher than those of the DR calculation, as 
in the constrained MSSM scenarios. Around the maxima of the curve, however, the difference 
between the two computations reaches 4-5 GeV. This indicates that for large mixing in the 
stop sector, X t ~ ^/6M$, the corrections which are formally of higher order (such as the 
three-loop terms controlled by the top Yukawa coupling and the strong coupling) are larger. 

We can now discuss the numerical results that we obtain with some of the spectrum 
generators in the three pMSSM scenarios defined in section I3~T1 Table HP shows the results 
obtained with the spectrum generators SPheno and SuSpect for the masses of the lighter and 
heavier CP-even Higgs bosons in the three pMSSM scenarios (the option for the pMSSM 
mode is not yet completely implemented in SoftSusy). It can be seen from the table that 
the results of the two codes agree within less than a GeV for Mh, as was the case in the 
constrained MSSM scenarios (see table [l}. The residual discrepancies are due, as discussed 
in section to the different ways in which the two codes implement the various threshold 
corrections to gauge and Yukawa couplings at the scale Q = M z 2 . On the other hand, it 
can be seen that the two codes are in excellent agreement for the heavier CP-even Higgs 
boson mass. This is due to the strict correlation between Mh and Ma'- contrary to the case 
of the constrained MSSM scenarios, where Ma is determined through the EWSB conditions 
and is thus sensitive to all of the details of the RG procedure, in the pMSSM scenario Ma 
is given as an input and is thus the same for both codes. 

The impact of the radiative corrections on the Higgs boson mass calculations is similar 
to what has been discussed in the case of the constrained scenarios and will not be repeated 
here. Furthermore, for the heavier CP-even Higgs boson, all of the theoretical uncertain- 
ties discussed in section 3.3 are very small, generating shifts of less than 100 MeV. In the 
following, we will therefore only concentrate on the uncertainties in the evaluation of the 
lighter Higgs boson mass Mh due to the renormalisation scheme, the renormalisation scale 
and the dependence on the external momentum. The results, obtained using SuSpect, are 
summarised in tabled] for the three pMSSM points. 



2 A further difference between the codes is that, in the pMSSM mode, SuSpect uses the physical Higgs 
boson masses in the computation of the various threshold corrections at Q = Mz , whereas SPheno switches 
to the physical masses only if one of the squared running masses becomes negative at Q = Mz- 
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Mass 


Code 


pMSSMl 


pMSSM2 


pMSSM3 


M h 


SPheno 


114.3 


118.8 


130.0 




SuSpect 


113.8 


118.4 


129.4 


M H 


SPheno 


1000.2 


1000.1 


999.4 




SuSpect 


1000.2 


1000.2 


999.5 



Table 11: The lighter and heavier CP-even Higgs boson masses, Mh and Mh, in the pMSSM 
scenarios as computed by SPheno2 .2.1 and SuSpect2 . 3. The SM and MSSM input param- 
eters are chosen as in section 13.11 



Shift 


pMSSMl 


pMSSM2 


pMSSM3 


Default value 


113.8 


118.4 


129.4 


FeynHiggs shift 


+1.6 


+ 1.9 


+4.4 


M EWSB = 150 GeV 
M E wsb = 2 TeV 


+1.1 
-0.7 


+ 1.5 
-0.8 


+4.4 
-0.8 


p 2 = in 1-loop 
Estimated 2-loop effect 


+1.6 
+0.3 


+1.8 
+0.2 


+3.0 
+0.2 



Table 12: Effect of renormalisation scheme, renormalisation scale and estimate of the external 
momentum dependence in the determination of Mh for the three pMSSM scenarios (all the 
shifts are in GeV). 

The first line of the table, after the default values given by SuSpect, shows the shift in 
the value of Mh when it is calculated in the OS scheme by the program FeynHiggs, where 
the weak scale input parameters are taken from SuSpect. As mentioned previously, the 
difference between the two calculations of Mh amounts to a couple of GeV in the first two 
scenarios, but reaches 4.4 GeV in the large-mixing scenario pMSSM3. 

To estimate the renormalisation scale dependence of the Higgs boson masses in the 
pMSSM scenarios, we can vary the scale M EW sb at which the corrections are computed. 
By default, in the pMSSM scenarios this scale is identified with the scale at which the input 
parameters are given. In order to estimate the related uncertainty, we evolve the parameters 
of the three pMSSM scenarios from the input scale Q = 1 TeV to some new scale, and 
compute again the Higgs boson masses in terms of the new values of the parameters (there 
is an option in SuSpect allowing for this study). The resulting shifts in the lighter Higgs 
boson mass for two alternative values of the scale, M EW g B = 150 GeV and M EWSB = 2 TeV, 
are shown in the middle part of table El As can be seen, the shifts in Mh are usually of the 
order of 1 GeV, but again reach 4.4 GeV in the case of the pMSSM3 scenario. 
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Finally, we see from the lower part of table^Jthat the external momentum dependence of 
Mh at the two-loop level, as estimated from the one-loop momentum dependence calculated 
in SuSpect, is expected to generate a shift of 300 MeV or less, as was approximately the 
case in the cMSSM scenarios discussed in the previous sections. 

3.5 Summary of the theoretical uncertainties 

From the results presented in the previous sections, it is possible to attempt an estimate of 
the theoretical uncertainty that still affects the calculation of the lighter MSSM Higgs boson 
mass. Such an estimate is relevant because it affects the information on the MSSM parameter 
space that can be extracted from the comparison between the theoretical prediction on Mh 
and the present experimental lower bounds. The phenomenological implications of this 
theoretical uncertainty will be further discussed in section 0] 

It has been shown in sections 13. 3. II and 3.4 that the discrepancies between the OS and DR 
calculations of the lighter Higgs boson mass Mh, which formally amount to a combination of 
two-loop effects controlled by the electroweak gauge couplings and three-loop effects, are of 
the order of 2-3 GeV in most of the considered cMSSM and pMSSM scenarios, but they can 
reach 4-5 GeV in the pMSSM scenario with large stop mixing. Another measure of the size 
of the uncomputed corrections is given by the residual renormalisation scale dependence, 
and in sections 13.3.21 and 13.41 it has been shown that the results of the DR calculation of 
Mh vary by 2-3 GeV when the renormalisation scale is varied in a reasonable range. The 
discrepancies between the results of the different spectrum generators in the constrained 
MSSM scenarios discussed in section 13.2.11 correspond generally to three-loop effects (i.e., 
two-loop differences in the determination of parameters that enter the Higgs mass corrections 
at one loop), and in the case of Mh are contained within a half GeV. Finally, another source 
of theoretical uncertainty is the approximation of zero external momentum employed in the 
two-loop corrections. In sections 13.3.31 and 3.4 it has been shown that this uncertainty in 
Mh is less than half GeV in the considered scenarios. 

The way these different results should be combined into a figure for the residual theoretical 
uncertainty on Mh is debatable, but we think that a reasonably conservative estimate of the 
global uncertainty can be taken as 3-5 GeV. The lower value is expected in most of the 
cMSSM and pMSSM scenarios, while the upper value is expected in some special scenarios 
such as in the pMSSM with large mixing in the stop sector. 

In Ref. |33] the combined effect of the corrections yet to be computed was estimated to 
induce an uncertainty in Mh of the order of 3 GeV. We find that our results are in good 
agreement with that earlier estimate for most of the scenarios that we have discussed in this 
paper. A recent analysis with the program FeynHiggs [IE], gives for the uncertainty on Mh 
a range that is similar to the one that we quote here, i.e. AM/, ~ 3-5 GeV, depending on 
the considered MSSM scenario. 

The fact that the predicted experimental accuracy in the determination of Mh at a future 
linear collider is 50 MeV [1H] indicates that a huge effort is still necessary in order to improve 
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the theoretical uncertainty on the Higgs boson masses up to the level required to compare 
with the forthcoming experimental results. An experimental effort is also needed to measure 
the SM input parameters more precisely, since they generate an error on Mh which is of the 
same order as the theoretical uncertainty, as will be discussed now. 

3.6 The impact of the experimental uncertainties 

It is important to discuss how the predictions for the MSSM Higgs boson masses are affected 
by the experimental errors in the determination of the input SM parameters. This is a 
source of uncertainty independent of the ones discussed above, and it can be reduced by 
more precise experimental information. 

It is well known, for example, that the precise values of the Higgs boson masses in the 
MSSM depend strongly on the exact value of the input top quark mass. Indeed, because the 
dominant one-loop corrections to the Higgs masses grow as M 4 4 , a shift of several GeV on 
M t will produce a significant variation of Mh- Since the experimental error on the top quark 
pole mass, AM t = ±4.3 GeV, is rather large, it will have a large impact on the numerical 
values of the Higgs masses. The benefits arising from a more precise measurement of M t for 
the determination of the MSSM parameters have been discussed in Ref. |50| . 

Tables IT3~1 and [T41 show the variation of Mh and Mh, respectively, in the six SPS scenarios, 
resulting from a ±ler variation in the top quark mass M t and in other SM input parameters, 
i.e. the bottom quark mass m(,(m(,) MS , and the strong and electromagnetic coupling constants 
a s (M z ) MS and a~^(Mz) MS , which are given in eqs. (JHJHQ; the other SM input parameters, 
such as Gf,Mz and M T , are so precisely measured that the impact of their experimental 
errors on the Higgs boson masses is negligible. 

As expected, the effect of the error on the top quark mass is rather large. In the case 
of the lighter CP-even Higgs boson, a 4.3 GeV variation in M t results in a shift in Mh of 
approximately 2 GeV in most SPS scenarios, 3 with larger (lower) values of M t increasing 
(decreasing) the obtained value for M h ; an exception is given by the scenario SPS5, where, 
due to a large stop mixing and a low tan (3 value, the top quark mass plays a more prominent 
role, and the shift can reach 3 GeV. In the case of the heavier CP-even Higgs boson, the 
variations are at the level of a few percent, except for the SPS4 scenario, where they reach 
the level of 10%. Differently from the case of the lighter Higgs boson mass, Mh can either 
decrease or increase with increasing (or decreasing) M t . This is due to the fact that, as 
discussed previously, the heavier Higgs boson masses are mostly sensitive to the values of 
the soft SUSY-breaking Higgs mass parameters, whose RGE evolution has a complicated 
dependence on the top Yukawa coupling. 

The impact on the value of Mh of the experimental error on the bottom quark mass, 
mb(mb) MS = 4.25 ±0.25 GeV, is significant only in the SPS4 scenario with high tan/5, where 

3 There is a well known rule of thumb for the so-called M™ aK scenario, according to which a variation of 
one GeV in M t results into a variation of one GeV in Mh- We notice that this rule does not generally hold 
in the constrained MSSM scenarios. 
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it can lead to a shift of about 100 MeV; otherwise, it is negligible. The effect, however, can 
reach the level of 10% or so for the determination of Mh in this high tan/3 scenario, since, 
in this case, the bottom Yukawa coupling plays an important role in the evolution of the 
soft SUSY-breaking Higgs parameters upon which Mjj depends. Note that lower (higher) 
values of the bottom quark mass in general increase (decrease) the values of Mh- 

The main effect of the error on the strong coupling constant is to generate an uncertainty 
in the top and bottom quark running masses, which propagates on the soft SUSY-breaking 
parameters and eventually on the Higgs mass spectrum. Varying a s (Mz) MS = 0.1172±0.002 
within its la range leads in general to an uncertainty of less than one per mille for Mh and 
one percent for Mh (an exception being the scenario SPS4 where the uncertainty for Mh is 
two percent). 

Finally, as one might have expected, the uncertainty arising from the error on the elec- 
tromagnetic fine structure constant, a _1 (M^) MS = 127.934 ± 0.027 (due to the hadronic 
contribution uncertainties), is essentially negligible: in all cases, it barely reaches the per 
mille level for both Mh and Mh- 

We now turn to the case of the pMSSM. First of all, in the considered scenarios the errors 
due to the variation of all SM parameters within their ler range are completely negligible 
in the case of the heavier CP-even Higgs boson: the maximal shift is less than 100 MeV 
(i.e. < 0.01%) and occurs when the top mass is varied within its 4.3 GeV uncertainty. As 
discussed previously, this is due to the fact that in the pMSSM the value M H is essentially 
controlled by the one of Ma, which is fixed to 1 TeV in these scenarios. Second, even in 
the case of the lighter Higgs boson, the errors due to the variation of a and m& are very 
small and can be safely neglected (the insensitivity to the precise value of m& is due to the 
moderate value chosen for tan/3 in the considered scenarios). 

The only SM parameters which have an impact on the value of Mh in the considered 
scenarios are M t and a s . The shifts due to their experimental errors in the three scenarios 
are shown in table ED The main effect of the variation of a s is to shift the DR top mass m t 
which enters the loop corrections to the lighter Higgs boson mass; the shift in Mh increases 
with the magnitude of the stop mixing and reaches the level of 300 MeV in the large mixing 
scenario. The largest uncertainty is, as expected, due to the experimental error on the top 
quark mass. The effect follows the same trend as in the SPS scenarios: for small mixing 
(pMSSMl) the induced error is of about 2 GeV, and it reaches the level of 4 GeV in the large 
mixing scenario (pMSSM3). To obtain an accuracy of 50 MeV on the lighter Higgs boson 
mass (as will be measured experimentally at the future colliders) one needs an experimental 
determination of M t to better than 100 MeV. Our results are in full agreement with Ref. [5D] . 
where the uncertainty in M h induced by the errors on M t , a s , m b and M w is discussed. 

In conclusion, the experimental uncertainties on the SM parameters lead to an error of 
about 3-4 GeV on the lighter Higgs boson mass Mh. This error is thus of the same magnitude 
as the total theoretical uncertainty which has been discussed in the previous section, but 
will be decreased when new measurements of the top mass become available. 
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Approximat ion 


SPSla 


SPS2 


SPS4 


SPS5 


SPS8 


SPS9 


Default 


112.1 


116.8 


114.1 


116.2 


115.4 


117.5 


AM t = ±4.3 GeV 
Am b = ±0.25 GeV 
Aq s = ±0.002 
Aa- 1 = ±0.027 


+1.59 
-1.55 

<0.01 

+0.02 

+0.09 
-0.06 

+0.02 
<0.01 


+2.12 
-2.04 

<0.01 
<0.01 

-0.06 

+0.07 

+0.01 
<0.01 


+1.73 
-1.69 

-0.07 
+0.09 

+0.06 
-0.05 

+0.01 
<0.01 


+2.09 
-3.05 

<0.01 
<0.01 

-0.09 
+0.03 

-0.01 
<0.01 


+2.26 
-2.14 

<0.01 
<0.01 

-0.05 
+0.05 

<0.01 
<0.01 


+ 1.89 
-1.84 

-0.01 
+0.01 

+0.09 
-0.10 

-0.01 
<0.01 



Table 13: The lighter CP-even Higgs boson mass, M^, in the six SPS scenarios, as computed 
by SuSpect using different input values for the SM parameters. The first row contains the 
default values in GeV, the other rows contain the shifts (in GeV) due to the different inputs. 



Approximation 


SPSla 


SPS2 


SPS4 


SPS5 


SPS8 


SPS9 


Default 


405.9 


1551.2 


354.4 


686.6 


550.3 


1056.6 


AM t = ±4.3 GeV 
Am b = ±0.25 GeV 
Aa s = ±0.002 
Aa- 1 = ±0.027 


+6.5 
-6.7 

-1.2 

+1.2 

+2.2 
-1.9 

+0.3 
-0.2 


+42.6 
-42.5 

-3.0 

+2.8 

-7.1 
+7.4 

+0.4 
-0.2 


+27.2 
-29.4 

-39.5 

+34.2 

+6.3 
-6.2 

+0.2 
-0.1 


-n.i 

+7.5 

-0.8 
+0.7 

+5.7 
-6.1 

+0.2 
-0.2 


+16.0 
-15.8 

-2.5 
+2.3 

+3.4 
-3.5 

-0.1 
+0.1 


-15.6 
+10.2 

-5.3 
+4.9 

+ 13.9 
-14.4 

-0.1 
<0.1 



Table 14: Same as table [13] for the heavier CP-even Higgs boson mass, M H . 



Approximation 


pMSSMl 


P MSSM2 


pMSSM3 


Default 


113.83 


118.39 


129.35 


AM t = ±4.3 GeV 
Aa s = ±0.002 


+2.31 
-2.21 

TO. 17 


+2.72 
-2.60 

T0.21 


+3.83 
-3.66 

T0.30 



Table 15: The lighter CP-even Higgs boson mass, Mh, in the three pMSSM scenarios, as 
computed by SuSpect using different input values for the SM parameters. The first row 
contains the default values in GeV, the other rows contain the shifts (in GeV) due to the 
different M t and a s inputs; the errors due to Am b and Aa are negligible. 



31 



Another independent source of uncertainty are the experimental errors in the determina- 
tion of the masses and couplings of the supersymmetric particles. A detailed study of this 
issue goes beyond the scope of our paper. 

4 Phenomenological consequences 

In this section, we discuss the phenomenological impact of the radiative corrections in the 
MSSM Higgs sector and of the related theoretical and experimental uncertainties. In par- 
ticular, we study the maximal values for the lighter Higgs mass as well as the constraints on 
tan/3 coming from the negative searches of the Higgs bosons at LEP2. The analyses will be 
performed in the context of the constrained MSSM models (mSUGRA, AMSB and GMSB), 
but also in the context of the unconstrained model (pMSSM). 

4.1 The upper bound on the lighter Higgs boson mass 

At the tree-level, the mass of the lighter MSSM Higgs boson Mh is bounded from above 
by Mz- Loop corrections increase this bound to ~ 135 GeV as discussed previously. The 
knowledge of the precise value of this upper bound is crucial for several reasons. 

- It is entirely conceivable that the lighter Higgs boson be discovered before the SUSY 
particles. In this case, it is important to accurately know the maximal possible value 
of Mh in order to discriminate between the SM and the MSSM. Furthermore, even 
if SUSY particles have been observed, the knowledge of this value, in the absence of 
extra information, could allow one to distinguish between the MSSM and some non 
minimal SUSY models where the Higgs sector is extended. 

- Even in the context of the MSSM itself, the knowledge of the maximal Mh value could 
allow one to discriminate between various scenarios of SUSY breaking. This approach 
would be then complementary to the study of the differences in the SUSY particle 
spectra induced by the various breaking schemes. 

The value Mh ~ 135 GeV has a rather special status in the decoupling regime Ma 3> 
Mz (which occurs very often in the constrained scenarios) where all of the MSSM 
Higgs bosons except for the lighter are too heavy to be observed experimentally at the 
next generation of high-energy colliders. Indeed, if Mh is below this critical value, the 
Higgs particle will dominantly decay into b quark pairs, while above this value it will 
dominantly decay into W bosons (for a detailed discussion of the decay modes, see 
Ref. [ST] for instance). 

4.1.1 The upper bound on M h in the pMSSM 

As discussed in section 3.4, in the pMSSM only a small subset of input parameters have a 
significant impact on the MSSM Higgs sector. In first approximation, to obtain the maxi- 
mal value of lighter Higgs boson mass one does not need to scan over the 22 free pMSSM 
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parameters, but, rather, to choose the relevant parameters in such a way that the one-loop 
radiative correction e in eq. (|llj) is maximised. In particular, one can obtain a reasonable 
approximation of the maximal Mh when one has: 

i) large values of the parameter tan (3, tan j3 > 20 (but still tan /3 < 60 to keep the bottom 
Yukawa coupling in a perturbative regime); 

ii) a decoupling regime with a heavy pseudoscalar Higgs boson, Ma ~ C(TeV); 

iii) heavy stops, i.e. large Ms values; we note, however, that heavier stops correspond to 
more fine tuning of the parameters in order to achieve the correct minimum of the 
Higgs potential [52] and we choose Ms = 2 TeV as a maximal value; 

iv) a stop trilinear coupling such that X t is close to +v / 6 M s , as exemplified by fig. El 

As a starting point, we realise the assumptions above by adopting the M™ ax scenario 
of Ref. jlHj, which was used as a benchmark point for the LEP2 Higgs analyses; we choose 
however to be conservative, scaling the relevant soft SUSY-breaking parameters by a factor 
of two and using the upper limit tan/3 ~ 60: 

x tan/3 = 60, M s = M A = 2TeV, A t = A b = V6M S , 

ftbench • M 2 ~ 2 M x = -/x = 400 GeV , M 3 = 0.8 M s . { ' 

For the central value M t = 178 GeV, and varying the values of a s and m& (recall that 
we are in the large tan/3 regime where the bottom Yukawa coupling is enhanced) in their ler 
allowed range, one obtains the maximal value of the lighter MSSM Higgs boson, 

M Cn X ch - 138 GeV for M t = 178 GeV. (14) 

This value is higher than the one which is often quoted in the literature, Mh ^ 135 GeV. 
The tendency of our DR calculation to give lower values of Mh than the OS calculation 
(using the same parameters as in Ref. jlS] we would find M h < 130 GeV) is compensated by 
the increase of the central value of the top quark mass and by the fact that we conservatively 
set the SUSY-breaking scale to 2 TeV. 

However, the bound in eq. (|14jl on the lighter Higgs boson mass is not yet fully optimised. 
In order to find the maximum maximorum 4 Mh value, one can still vary in a reasonable 
range the SUSY parameters entering the radiative corrections and add the estimated theo- 
retical (and experimental) uncertainties. We have therefore studied how to tune the pMSSM 
parameters in such a way that the corrections to Mh are maximised. 

For instance, the one-loop electroweak corrections involving Higgs, gauge bosons and 
their fermionic superpartners tend to lower the value of Mh when the chargino and neutralino 
masses are large. On the other hand, those corrections get small, allowing for a maximal Mh, 
when the chargino and neutralino masses are of the same order of the gauge boson masses 
(see e.g. Ref. [TH] for a discussion). For fi ~ M 2 ~ |Mi ~ 150 GeV, the upper bound on M h 

4 We thank Daniel Treille for suggesting this expression. 
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can be pushed upwards by about 2 GeV, while the lightest chargino mass is still larger than 
104 GeV, as required by the negative searches at the LEP2 [46J. 

The gluino mass enters the two-loop QCD corrections in the Higgs sector and can have 
a non-negligible impact on the value of Mh- For large tan/3, lowering the value of the gluino 
mass parameter to M 3 ~ 0.4 Mg leads to an increase of Mh of about 2 GeV. 

Another increase of the upper bound on Mh can be obtained by relaxing the equality of 
the stop mass parameters at the weak scale (this turns out to affect Mh through changing 
the threshold correction to the top Yukawa coupling). For large tan/?, assuming Mj ~ 
2Mf ~ 2.5 TeV (which still keeps the EWSB scale, defined as the geometric mean of the 
two stop masses, in the vicinity of 2 TeV), one obtains an increase of Mh of about 2 GeV. 

One can then optimise the soft SUSY-breaking parameters to maximise all the effects 
mentioned above. Using the codes SPheno and SuSpect, we have performed a scan of 
the MSSM parameter space with the following constraints: (i) Mewsb < 2 TeV, even if 
some SUSY particles (in particular the heavier stop) may have masses above this value; 
(ii) the trilinear couplings of the third generation sfermions should not exceed the value 
At,b,r ~ 3Mewsb, from which problems with charge and colour breaking (CCB) minima 
might appear [HHj : (iii) all sparticles, in particular the lightest chargino and the gluino, 
should have masses that exceed their experimental bounds; (iv ) the resulting spectrum 
should lead to acceptable contributions to the p parameter, to the anomalous magnetic 
moment of the muon (g — 2) M and to the radiative b — > s'y decay 5 . The two codes give rather 
close results, the differences being of the order of 1 GeV. To provide a conservative upper 
bound on Mh, we quote the result of the code that gives the largest output value. 

Using the central value of the top quark mass, we find the following upper bound on the 
lighter Higgs boson mass in the unconstrained MSSM: 

M h < 143 GeV for M t = 178 GeV. (15) 

Thus, tuning the SUSY parameters 6 in the way discussed above is very important, leading 
to an increase of ~ 5 GeV of the upper bound on Mh- 

If we vary the top quark mass within its la experimental range we obtain the following 
upper bounds on Mh for the lower and upper M t values 

M h < 138 GeV for M t = 173.7 GeV, 

M h < 148 GeV for M t = 182.3 GeV. (16) 

Finally, to obtain the most conservative maximal value of Mh, one has to include the 
theoretical uncertainty. For the highest value of the top quark mass, M t = 182.3 GeV, when 

5 The (g — 2) M bound is easily satisfied as the second-generation sfermion masses can be taken to be large 
without affecting the Higgs mass calculation, while the Ap constraint is in general weaker than the related 
CCB constraint. On the other hand, constraints coming from dark matter are not studied, being beyond 
the scope of this paper. 

6 For files detailing the choices of parameters that lead to the maximal Mh and the resulting MSSM mass 
spectra, see the SuSpect jS] and SPheno 10, webpages. 
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we also include the theoretical uncertainty that we estimate to be AM/, « 4 GeV in this 
case, we obtain the following maximum maximorum bound on the lighter Higgs boson mass 
in the unconstrained MSSM: 

M h < 152 GeV. (17) 

Results for the maximal Mh value in the unconstrained pMSSM, based on the two-loop 
OS calculation implemented in FeynHiggs, have been presented e.g. in Ref. [33], where the 
upper bound M h < 140 GeV can be inferred for M s = 2 TeV and M t = 179.4 GeV. A 
recent update [3BJ, using the latest upper value M t = 182.3 GeV, gives Mh ^ 143 GeV. 
In both cases, the MSSM parameters are chosen as in the M™ 5 ™ scenario of Ref. |45j . and 
the theoretical uncertainty is not taken into account. Including the theoretical uncertainty, 
these results become only a few GeV lower than the upper bounds discussed above. 

4.1.2 The upper bound on Mh in the constrained MSSM 

In the mSUGRA, GMSB and AMSB scenarios, the various parameters which enter the 
radiative corrections are not all independent, due to the relations between SUSY breaking 
parameters that are set at the high-energy scale. In addition, the radiative EWSB constraint 
must be fulfilled for each set of input parameters (in the pMSSM, this is automatic since 
Ma and ii are used as input parameters). Thus, in contrast to what occurs in the pMSSM, 
it is not possible to freely tune all relevant weak-scale parameters in order to get a maximal 
value of Mh- This should in principle make the analysis more complicated. However, since 
we have only a small set of input parameters in the constrained case, one can easily perform 
scans of the parameter space. 

For the purpose of illustration, we show in figs. EHH1 how Mh depends on some relevant 
parameters of the mSUGRA, GMSB and AMSB scenarios. The results shown in the figures 
are obtained with Sof tSusy. Fig. 0] shows the value of Mh as the background density in 
the plane Aq — tan/3 of the mSUGRA model, itlq and m 1( / 2 are set to 1 TeV in order to 
have an average stop mass Ms of the order of 2 TeV, and \x is taken positive. In the white 
region in the upper-left corner of the plot a sparticle squared mass is negative, resulting 
in an unacceptable scalar minimum. Contours of M$ and Mh are also shown. Large stop 
mixing is achieved for large and negative A , around —4 TeV, corresponding to the maximal 
Mh ~ 128.5 GeV region in the figure. 

Fig. El shows the value of Mh as the background density in the plane A — tan (3 of the 
GMSB model. The messenger mass M mess is taken equal to 10 5 x A. Moreover, N meas = 8 
and [a is taken positive. It can be seen that the maximal value for Mh, which is around 
123.5 GeV in this example, is obtained for large values of the SUSY-breaking scale A, which 
correspond to large stop masses. 

Finally, in fig. E]we show Mh in the 7713/2 — tan/5 plane of the AMSB model. The gravitino 
mass m 3 /2 is varied in a range that allows for stop masses of the order of 2 TeV, the scalar 
mass term m is set to 1 TeV, and ii is taken to be negative. Even in this case, it can be 
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Figure 4: The lighter Higgs boson mass M^ in the tan f3 — A plane of the mSUGRA scenario. 
Other input parameters are m = mi/2 = 1 TeV and /i > 0. The SM input parameters are 
fixed to their default values. The lighter Higgs boson mass is shown as background density, 
as measured by the bar on the right. Dashed contours are lines of equal Ms, dotted contours 
are lines of equal Mh- The white region contains scalars with negative squared mass. 

seen that Mh increases for larger values of m 3 / 2 , which correspond to larger values of the 
stop masses. 

Plots like the ones discussed above can give useful indications about the choice of param- 
eters that leads to the maximal Mh- However, in order to get a reliable determination of the 
maximal Mh in a given SUSY-breaking scenario it is necessary to scan through the allowed 
range of values for all of the relevant parameters. Such scans were previously performed, e.g., 
in Refs. |3BJ EH EH] • The most recent among these analyses, Ref . J3H] , used the combination 
of FeynHiggs with different private codes jSl] for the computation of the MSSM parameters 
in the various SUSY-breaking scenarios. For convenience, we take the same ranges of soft 
SUSY-breaking input parameters as in Ref. [3%] : 

mSUGRA: 50 GeV < m < 1 TeV, 50 GeV < m 1/2 < 1 TeV, \A \ < 3 TeV; 
GMSB: 10 TeV < A < 200 TeV, 1.01 < M mcss /A < 10 5 , 1 < V mcss < 8; 



AMSB: 20 TeV < m 3/2 < 100 TeV, 50 GeV < m < 2 TeV. 

Moreover, we require 

1 < tan/3 < 60 

and we allow for both signs of \l. The SM input parameters are fixed to their default values 
as in eqs. (JBJ)-©- In particular, we take M t = 178.0 ± 4.3 GeV |2S], whereas Ref. [SH] used 
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Figure 5: Same as fig.HJin the tan/3— A plane of the GMSB scenario. Other input parameters 
are M mess = 10 5 x A, A" mess = 8 and fi > 0. 



the old world average M t = 175 GeV. 

To avoid the need for excessive fine-tuning in the EWSB conditions, we also impose an 
additional bound on the weak-scale parameters, i.e. that the geometrical average of the stop 
masses is less than 2 TeV, 



Ms 



Mewsb = y/ r m il m i 2 < 2 TeV • 



However, it is important to bear in mind that, in the absence of a compelling criterion 
to define the maximal acceptable amount of fine-tuning, the choice of the upper bound on 
Ms is somehow subjective. This upper bound on Ms affects the maximal attainable value 
of M h , as larger values of M s push the value of M h upward. 

Taking the range of SUSY-breaking parameters described above and scanning with 
SoftSusy, SPheno and SuSpect leads to the results for the maximal M^ given in table 
IT6| where M t is set to 173.7, 178.0 and 182.3 GeV; we also display, for completeness, the 
upper bounds in the case of the pMSSM. 

We display in table HB] only the result of the code that gives larger upper bound on M^, 
the differences between the results of the three codes in the search for the maximal M^ can 
amount to 1 GeV, somewhat larger than those found in section 13.2. II for the SPS scenarios. 
Apparently, the higher order effects that induce the residual differences among the codes 
are enhanced for the extreme choices of the SUSY parameters that give rise to the maximal 
values of M^. In addition, in the last column of the table we show our most conservative 
estimate for the upper bound on the lighter Higgs boson mass in the mSUGRA, GMSB 
and AMSB scenarios. For this, we take M t = 182.3 GeV, the central values for the other 
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Figure 6: Same as fig. 0] in the tan/3 — m plane of the AMSB scenario. Other input 
parameters are 777.3/2 = 100 TeV and \i < . 

SM input parameters and add to the results an uncertainty of approximately 4 GeV (which 
includes both the theoretical uncertainty as estimated previously, and the much smaller error 
due to the variation of the remaining SM input parameters). 

It can be seen that, for the central value of M t , the GMSB model has the tightest bounds 
upon the maximal lighter Higgs mass, M^ < 123.7 GeV, followed closely by the AMSB 
model, where the maximal value is obtained for jjl < (for the central value of the top mass, 
setting /i > would induce a 0.6 GeV decrease in the maximal Mh). The mSUGRA model 
has a much looser bound, ~ 5 GeV above the GMSB case, but this bound is more than 10 
GeV smaller than in the case of the pMSSM. This pattern can be qualitatively understood by 
considering in each model the allowed weak scale values of A t , which essentially determines 
the stop mixing parameter X t and, therefore, the value of Mh- The GMSB model has A t = 
at the relatively low scale M mess < 2 x 10 7 TeV, thus its magnitude does not increase greatly 
in the RG evolution down to Mewsb- A small value of A t implies a smaller value for the stop 
mixing, therefore producing less high Mh (see eq. (IlljO . The AMSB model has a non-zero 
At that is fully predicted at any renormalisation scale in terms of the Yukawa and gauge 
couplings, and of the overall SUSY breaking scale m 3 / 2 which determines also the common 
squark mass Ms- Thus, the ratio between X t and Ms cannot be tuned at will in order to get 
the maximal value of Mh- Nevertheless, X t /Ms turns out to be sizeable at the weak scale, 
resulting in a larger value of M^ than in the GMSB model. Finally, in the mSUGRA model 
the only restriction on A t (Mewsb) comes from \A \ < 3 TeV, allowing for close to maximal 
mixing, but only for values of the SUSY scale lower than the value Ms = 2 TeV which is 
used as input in the pMSSM case; if we allow for larger and negative values of A we can 
get closer to, but not saturate, the pMSSM bound. In fact, it can be seen from fig. El that, 



38 



M max 


M t = 173.7 GeV 
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= 178.0 GeV 
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= 182.3 GeV 


conservative 


pMSSM 


138 


143 


148 


152 


mSUGRA 


126.2 




129.0 




131.7 


136 


GMSB 


120.8 




123.7 




126.7 


131 


AMSB 


122.0 




124.6 




127.1 


131 



Table 16: Maximal lighter Higgs boson mass (in GeV) in the pMSSM, mSUGRA, GMSB 
and AMSB scenarios as obtained by the three spectrum generators for three values of the 
top quark mass. The last column gives the most conservative estimate, where M t = 182.3 
GeV and a theoretical uncertainty around 4 GeV is added. 



in the unconstrained MSSM, the maximal Mh is obtained for positive values of X t . On the 
other hand, in the constrained models the RG evolution tends to drive A t (thus, X t ) towards 
negative values. In order to get at the weak scale a positive X t large enough to obtain the 
maximal Mh, we should start from an unreasonably large positive value of Aq. 

Note that table El also shows that changing M t by one standard deviation adds or 
subtracts about 3 GeV to the maximal Mh value in each of the cMSSM scenarios, to be 
compared with the 5 GeV variation discussed previously for the pMSSM. 

In Ref jHHl, the upper bounds M h < 126.6,123.2 and 124.5 GeV are quoted for the 
mSUGRA, GMSB and AMSB models, respectively. These results follow the same qualitative 
pattern as those given in table El However, the good agreement found in the GMSB and 
AMSB scenarios using the central value of M t is to some extent accidental: as discussed 
above, our DR calculation gives results for Mh that are a few GeV smaller than the ones of 
the OS calculation in FeynHiggs, but this effect is compensated by the higher central value 
of M t taken here (M t = 178 GeV instead of M t = 175 GeV). Also, the authors of Ref. [Ml 
did not impose an upper bound on Ms, resulting for example in stop masses up to 10 TeV 
in GMSB (and therefore an enhancement of Mh). Finally, the case \i < was not considered 
in their analysis, which resulted in a reduced upper bound on Mh in the AMSB case. 

Note that the bounds on Mh discussed above will not be in conflict with any phenomeno- 
logical constraint. Since the obtained SUSY spectra in the cMSSMs is rather heavy, one 
easily evades the collider experiment bounds on the sparticle and Higgs boson masses and 
satisfies the indirect constraints from from high-precision measurements such as (g — 2) M , 
BR(6 — > 57) and Ap (the dark matter constraint is not included). 



4.2 The lower bound on tan/3 

The present LEP2 limit on the mass of the SM Higgs boson, M#o > 114.4 GeV jSHj, sets 
strong constraints on the parameter space of the MSSM Higgs sector. In general, a detailed 
investigation of the Higgs boson production and decay modes would be required to establish 
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whether a point in the MSSM parameter space is excluded by the experimental searches. 
However, in the decoupling regime where Ma 3> Mz, the lighter CP-even Higgs boson of the 
MSSM has the same couplings to fermions and gauge bosons as the SM Higgs boson (and 
hence the same production and decay properties) thus the LEP2 bound can be extrapolated 
to Mh > 114.4 GeV. The tree-level value of Mh decreases for small values of tan/3, as 
can be seen from eq. (|10J1. and the lower experimental bound on Mh translates into a lower 
bound for tan/3. 

In this section, we provide such a bound in the different SUSY breaking scenarios dis- 
cussed previously. Such investigation is interesting for several reasons: 

- An appealing possibility is that the Yukawa couplings of the third generation fermions 
are unified [e.g., hb(M GUT ) = h T (M GUT ) in minimal SU(5)]. In such models, the value 
of tan /3 is of critical importance and determines the viability of bottom-tau Yukawa 
coupling unification 



The experimental determination of the value of tan (3 is an important measurement to 
be performed once the supersymmetric and Higgs particles are detected. Unfortunately, 
such a measurement can be difficult (see Ref. JHH] for instance): in most cases the 
relations between the SUSY or Higgs parameters or observables which can be used 
for this determination depend on cos(2/3) rather than on tan/3 itself. This quantity 
becomes approximately equal to unity as soon as tan/3 > 3, leaving little room for 
an accurate measurement of tan/3, if low values of tan/3 are excluded by the LEP2 
constraint 8 . It is thus interesting to see whether there are still viably low values of 
tan/3 that allow for a measurement. 

- A related reason is that, in the MSSM Higgs sector, the parameter space of the decou- 
pling regime expands for higher values of tan /3. For tan /3 > 5 for instance, the lighter 
Higgs boson becomes SM-like as soon as Ma exceeds the maximal value of Mh- This is 
not the case for low values of tan/3: even for A boson masses of several hundred GeV, 
the lighter MSSM Higgs boson has different couplings to fermions and gauge bosons 
(and thus different decay and production properties) than the SM Higgs boson. 

4.2.1 The lower bound on tan/3 in the pMSSM 

In the pMSSM, the absolute lower limit on tan/3 can be obtained by maximising the lighter 
Higgs mass in such a way that it exceeds the value Mh > 114.4 GeV [55] from direct 
searches at LEP2 (for earlier studies, see e.g. Refs. [HS1ES1)- Therefore, one can perform the 
same analysis as for the determination of the maximal Mh value discussed in the previous 

7 Also, for tan/3 < 4, one is close to the infrared quasi fixed point regime of the MSSM |5S], where the 
lighter Higgs boson mass M/, becomes smaller j57j . 

8 Note that, for high values, tan/3 could be possibly measured by considering associated production of the 
heavier MSSM Higgs bosons with bottom quarks |Bjj], tau fusion into Higgs bosons [HI] or stau decays into 
Higgs bosons |rJ2] . 

40 



subsection, except that here, one looks for the minimal value of tan/5 which leads to Mh > 
114.4 GeV. Thus, one has to be in the regime defined by the items (ii)-(iv) of section 4.1, 
i.e. being in the decoupling regime with a large mixing scenario and using large stop masses. 
In addition, as in the case of the maximal Mh, one still has to tune the values of /x and of 
the soft SUSY-breaking parameters in order to maximise the radiative corrections. 

To obtain the lower bound on the parameter tan /3, we have performed a scan similar to 
the one discussed in section 4.1.1, looking for the maximal value of Mh corresponding to a 
given value of tan/?. In fig. [?l we display the variation of the maximal Mh in the pMSSM 
as a function of tan/3. The dotted, full and dashed lines show the values of Mh for the top 
mass values M t = 173.7, 178.0 and 182.3 GeV, respectively, while the dash-dotted line on 
the top is for the conservative case where M t = 182.3 GeV is used and a 4 GeV theoretical 
uncertainty is added to Mh. 
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Figure 7: Upper bound on the lighter Higgs boson mass Mh in the pMSSM as a function of 
tan/3, as obtained from a full scan of the parameter space. The dotted, full and dashed lines 
correspond to the top mass values M t = 173.7, 178.0 and 182.3 GeV, respectively, while the 
dash-dotted line on the top is for the conservative case where M t = 182.3 GeV is used and 
a 4 GeV theoretical uncertainty is added to Mh- 



As can be seen, for the default value M t = 178 GeV the LEP2 bound of 114.4 GeV on 
Mh is always satisfied and therefore, for this top quark mass value, no empirical bound on 
tan/3 (provided that it is larger than unity) can be derived in the pMSSM. This is of course 
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also the case for the larger value of M t = 182.3 GeV, and a fortiori for the conservative 
case where a theoretical uncertainty is taken into account, where all values 1 < tan/? < 60 
are allowed by the LEP2 constraint. Only in the case of a lighter top quark, M t = 173.7 
GeV, the range 1.1 < tan/3 < 1.6 is excluded if the LEP2 bound M h > 114.4 GeV is to 
be satisfied. However, if a theoretical uncertainty of 4 GeV on Mh is included (meaning, 
in practice, that the LEP2 Higgs mass bound translates to the bound Mh > 110.4 GeV 
on the prediction obtained without including the theoretical uncertainty) no bound on the 
parameter tan/3 can be obtained from the LEP2 constraint. 

The local maximum around tan/3 = 1.2 for the maximal Mh in fig. 0is due to the fact 
that the top Yukawa coupling is maximised by the threshold corrections in that point. On 
the other hand, in the large tan/3 regime one recovers the maximal Mh values discussed in 
the previous section, as shown in the right-hand side of the figure. In particular, one sees 
that the upper bound on the lighter Higgs boson, in the conservative case where M t = 182.3 
GeV and a 4 GeV theoretical uncertainty on M h is added, is slightly above 150 GeV. 

4.2.2 The lower bound on tan/3 in the constrained MSSM 

In the constrained MSSM scenarios, we parameterise the 95% C.L. bounds from LEP2 in 
the sin 2 (/3 — a)-Mh parameter space [55, , rather than just using the lower bound on M^, 
in case the decoupling regime is not obtained. In practice, the lowest values of tan/3 that 
saturate the LEP2 limit are indeed obtained in the decoupling regime. To reduce the volume 
of the space that is scanned over, we use the parameter ranges as described in section l4~D for 
the mSUGRA, GMSB and AMSB cases. The results are summarised for these three cases 
in table IT7I where the values M t = 173.7, 178.0 GeV and 182.3 GeV are used, while for the 
other SM inputs we used the default values given in eqs. ©-©• The last column shows 
the most conservative lower bounds, obtained by adding a theoretical uncertainty of 4 GeV 
to the prediction of the lighter Higgs boson mass. For completeness, we also include in the 
first row the results in the pMSSM case discussed previously. In each entry of the table, we 
quote the result obtained with the code that gives the weakest lower bound on tan/3. 



tan/3 min 


M t = 


= 173.7 GeV 


M t = 


= 178.0 GeV 


M t = 


= 182.3 GeV 


conservative 


pMSSM 


1.6 


- 


- 


- 


mSUGRA 




2.8 




2.4 




2.1 


1.9 


GMSB 




4.2 




3.3 




2.7 


2.2 


AMSB 




3.7 




3.1 




2.7 


2.3 



Table 17: Lower bounds on the parameter tan/3 in the pMSSM, mSUGRA, GMSB and 
AMSB scenarios for three values of the pole top mass. The last column gives the most 
conservative estimate, where we use M t = 182.3 GeV and add a 4 GeV theoretical uncertainty 
on M h . 
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One can see from the table that, as expected, the minimal values of tan/3 follow the same 
pattern as in the evaluation of the maximal Mh value. The bounds are most stringent in 
the GMSB scenario, closely followed by the AMSB scenario, and finally by the mSUGRA 
scenario. In all cases the bounds are much stronger than in the pMSSM. In fact, contrary 
to the pMSSM case, in the constrained scenarios there is always a significant bound on 
the parameter tan/? even for the largest value of M t , tan/3 > 2.1-2.7, and when a 4 GeV 
theoretical uncertainty on M h is included, tan/3 > 1.9-2.3. This is merely a consequence 
of the ~ 14 GeV difference between the maximal values of Mh in the constrained and 
unconstrained cases. 

In Ref. [SB! the lower bounds tan/3 > 2.9, 3.1 and 3.8 for the mSUGRA, GMSB and 
AMSB scenarios, respectively, are quoted. In contrast to our findings, summarised in table 
IT7| the bound quoted in Ref. [HE] for the AMSB scenario is stronger than the one quoted 
for the GMSB scenario, presumably due to the fact that the authors did not consider the 
case ft < in their analysis 9 . On the other hand, our results for the mSUGRA and GMSB 
scenarios, using the central value of M t , are relatively close to the results of Ref. (SHJ- This 



is again due to a compensation between the tendency of the DR calculation to give lower 
values of Mh than the OS calculation (thus, stronger bounds on tan/3), and the increased 
value of M t , resulting into higher values of M h (thus, weaker bounds on tan/3). 

4.3 Implications for the decays of the Higgs bosons 

As discussed in the previous subsections, in the most conservative case where the top quark 
mass is set to M t = 182.3 GeV and a theoretical uncertainty of 4 GeV is added to M h , the 
upper bound on the lighter Higgs boson mass in the unconstrained MSSM can reach the 
level of 150 GeV, and there is no lower bound on the parameter tan/3 from the negative 
Higgs boson searches at LEP2. Such large values of Mh and small values of tan/3 would have 
a decisive impact on the decays of the MSSM Higgs bosons (and also on their production 
rates) and thus on their experimental detection at future high-energy colliders, in particular 
at the LHC, as will be briefly discussed below. 

In the decoupling regime, the lighter Higgs boson h behaves approximately like the SM 
Higgs particle. In particular, it has the same couplings to fermions and gauge bosons, thus 
the same branching fractions and production cross sections, provided that the SUSY particles 
are heavy enough to make their contribution to the loop induced decays safely negligible [S5] . 
In fig. |Hl obtained with the code HDECAY, we display the branching ratios of the SM Higgs 
boson (equivalent to those of h in the decoupling regime) as a function of its mass. 

As can be seen, while the h —>■ bb decay mode is dominant for masses below Mh ~ 135 
GeV, it is the decay h — > WW (with one of the gauge boson being off-shell) which becomes 
the most important one above this Higgs mass value. In particular, for Higgs boson masses 



9 We have checked that in the considered scenarios the SUSY contribution is within the 2a range for 
(g — 2) p (see Ref. ]§%\ for a recent analysis). 
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close to the upper bound, Mh ~ 150 GeV, the decay into W boson pairs reaches a branching 
ratio of the order of 70%. Contrary to the bb decay mode, which suffers from a huge hadronic 
background at the LHC, the WW decay mode is much cleaner if one considers the leptonic 
decays of one or two W bosons, and therefore becomes the most promising detection channel 
at the LHC gEj. 

Because the partial decay width of the lighter Higgs boson into WW final states is 
enhanced by the gain in phase space, almost all of the other decay modes will be suppressed. 
The relatively clean decays into r-lepton pairs, which are a useful signal when the lighter 
Higgs boson is produced in the fusion of two gauge bosons at the LHC pp — > qqh with 
h — » t + t~ |66j, can be suppressed to the level of 2% and this detection channel might 
therefore be disfavoured. The rare but very clean decay of the lighter Higgs boson into two 
photons is not strongly affected, being at the level of a few per mille, and still constitutes a 
viable detection channel 10 . 

Another consequence of a large Mh would be that, because of the increase of the phase 
space, the branching ratio for the decay of the lighter Higgs boson into Z boson pairs can 
reach the level of 10%. This would make the detection channel h — > ZZ* — > 4^ ± more 
important at the LHC, and the larger number of events could eventually allow for the 
determination of the spin-parity quantum numbers of the Higgs particle from the study of 
the angular correlations between the final state leptons |57j . 
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Figure 8: The branching fractions of the SM Higgs boson (equivalent to those of the lighter 
MSSM Higgs boson h in the decoupling regime) as a function of its mass. 

Let us now turn to the consequences of the constraints on the parameter tan/5 obtained 
from the negative searches of the MSSM Higgs bosons at LEP2. First of all, we recall that the 



10 The branching ratio for this decay is not suppressed because it is dominantly mediated by W boson 
loops, whose contribution becomes maximal near the WW threshold. 
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lower tan (3 bounds discussed in the previous subsection, have been derived in the decoupling 
regime where the pseudoscalar Higgs boson is rather heavy, Ma = 2 TeV. To derive this 
constraint outside the decoupling regime, one can no longer use the LEP2 bound on the SM 
Higgs boson mass, but rather, the sin 2 (/3 — a)-Mh parameter space excluded by the LEP 
collaborations [HE] in the search of the MSSM Higgs bosons in all production channels. A 
detailed analysis of the lower bound on tan/3 for any value of Ma and its impact on the 
decays of the MSSM Higgs bosons and on their production rates at the various high-energy 
colliders is beyond the scope of this paper. Nevertheless, a rough analysis shows that, in the 
conservative case discussed above, values of tan/3 ~ 1.8 are still allowed for Ma values as 
low as 200 GeV, which leads to the following qualitative conclusions: 

- The decays of the three neutral MSSM Higgs bosons into bb and t + t~ final states may 
not be too strongly enhanced. The lighter MSSM Higgs boson will then behave more 
like the SM Higgs particle even outside the decoupling regime. 

The decays of the heavier CP-even Higgs bosons into two gauge bosons and into two 
lighter Higgs bosons would be less suppressed. In particular, if the decay H — > hh is 
allowed it could lead to the possibility of measuring the trilinear H hh coupling. 

The decays of the pseudoscalar and charged Higgs bosons into a gauge boson and the 
lighter Higgs boson, A —>■ hZ and H ± — > hW^, could provide interesting detection 
signals. 

A detailed analysis of these features will be presented in a forthcoming publication 



5 Conclusions 

In this paper, we provided a detailed analysis of the Higgs sector of the Minimal Supersym- 
metric extension of the Standard Model, both in the constrained cases such as mSUGRA, 
GMSB and AMSB and in the unconstrained case, where 22 parameters are allowed to be 
free at the EWSB scale. We have provided a very precise determination of the masses of the 
neutral Higgs bosons including all the presently available radiative corrections and using the 
latest value of the top quark mass, as recently measured at the Tevatron. 

In three public RGE codes which evaluate the supersymmetric particle spectrum, Sof tSusy, 
SPheno and SuSpect, we have implemented all of the available radiative corrections for the 
Higgs boson masses and for the effective potential, in the DR renormalisation scheme which 
is more natural for these spectrum generators. In particular, we have incorporated the full 
set of one-loop radiative corrections to the Higgs self-energies and to the tadpoles, includ- 
ing the contributions of all SUSY particles and taking into account the external momentum 
dependence. We have also included the dominant two-loop radiative corrections at vanish- 
ing external momentum, i.e. the corrections which involve the strong gauge coupling and 
the Yukawa couplings of the third generation fermions. To complete this picture, we have 
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derived the missing radiative corrections controlled by the r lepton Yukawa coupling, which 
however turned out to be numerically small. 

We have performed a detailed comparison of the three programs, in six SPS scenarios for 
the mSUGRA, GMSB and AMSB cases as well as in three scenarios for the unconstrained 
MSSM. We have found a very good agreement among the results of the three codes: to 
within half GeV for the mass of the lighter CP-even Higgs boson and within 1% for the 
heavier CP-even and CP-odd Higgs bosons. We therefore conclude that the three codes 
provide now a reliable and very precise determination of the neutral Higgs boson masses, 
both in constrained and unconstrained MSSM scenarios. 

We have then studied the effects of the various radiative corrections on the Higgs boson 
masses in constrained scenarios and in the DR scheme. The corrections due to the strong 
interactions, 0(a t a s + a&a s ), increase the lighter Higgs boson mass by up to 6 GeV, while 
the pure Yukawa corrections due to the top and bottom quarks, 0(af + a t ai> + af), can 
decrease it by 1 GeV; the corrections due to the r- Yukawa couplings, 0{a^ + af,« T ), can 
be safely neglected. The impact of some higher order effects, such as those originating from 
the inclusion of the two-loop standard QCD and one-loop electroweak corrections to the 
top quark mass, the resummation of the SUSY-QCD corrections in the bottom quark mass, 
the use of the two-loop RGEs for the sfermion masses and the use of the pole Higgs boson 
masses in the one-loop corrections, has also been analysed. These higher order effects can 
alter the results by up to half a GeV for Mh and up to a few percent in Mh- 

A detailed study of the theoretical and experimental uncertainties in the determination 
of the CP-even Higgs boson masses has been performed in both the constrained and uncon- 
strained MSSM scenarios. The comparison between the three RGE codes and the program 
FeynHiggs, where the MSSM Higgs boson masses are evaluated in the On-Shell renormalisa- 
tion scheme, gives an estimate of the scheme dependence of the predictions. This dependence 
is in general of the order of 2 GeV, but can reach the level of 4-5 GeV in some cases, such 
as in the pMSSM scenario with large mixing in the stop sector. The dependence on the 
scale at which the effective scalar potential is minimised and the radiative corrections are 
computed is significantly improved by the inclusion of the two-loop radiative corrections, 
leaving a variation of less than 3 GeV on Mh (apart from the case of the pMSSM scenario 
with large stop mixing, where the variation is larger). The effect of approximating the ex- 
ternal momentum in the two-loop Higgs boson self-energies to zero rather than to the mass 
of the Higgs boson has been estimated to lead to a shift of less than half a GeV on Mh and 
to an even smaller shift on Mh- 

All of these shifts indicate the size of the higher order corrections that are left uncom- 
puted, thus they are a measure of the theoretical uncertainty which affects the calculation. 
As a summary, we conclude that these effects induce a theoretical uncertainty on the mass 
of the lighter Higgs boson of 3-4 GeV in the constrained MSSM and up to 5 GeV in the 
unconstrained MSSM, in particular when the mixing in the stop sector is such that Mh is 
maximised. The uncertainties on Mh due to the experimental errors in the measurement 
of the SM input parameters are approximately of the same size. Indeed, a variation of the 
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top quark mass within its la range, M t = 178.0 ± 4.3 GeV, leads to a shift of about ±2 
to ±4 GeV on Mh, depending on the scenario. The shifts on Mh from the errors in the 
measurement of mb and a s are less than half a GeV, while the error on a has a negligible 
impact. Thus, both a theoretical and an experimental effort are needed to attain the level 
where the lighter Higgs boson mass can be predicted with a precision of less than 100 MeV, 
i.e. the level of precision at which this mass can be measured at the future colliders. 

Finally, using the RGE codes, we have performed a phenomenological study of the impact 
of the radiative corrections, as well as the theoretical and experimental uncertainties, on the 
Higgs sector of the constrained and unconstrained MSSM. In particular, using the latest 
Tevatron value of the top quark pole mass, M t = 178.0 ± 4.3 GeV, we have derived the 
maximal value of Mh in the MSSM, as well as the minimal value of the parameter tan/5 that 
is allowed by the negative search of the MSSM Higgs bosons at LEP2. In the unconstrained 
MSSM, we get an upper bound Mh < 143 GeV and we show that all values of tan/3 are 
allowed by LEP2 constraints. In the constrained MSSM scenarios, as a result of the chosen 
ranges for the soft SUSY-breaking parameters and of the relations among them, the bounds 
on M h are stricter than in the pMSSM: 129.0, 123.7 and 124.6 GeV in the mSUGRA, 
GMSB and AMSB cases respectively. Similarly, in the constrained scenarios it is possible 
to derive lower bounds on tan/?: 2.4, 3.3 and 3.1 in the mSUGRA, GMSB and AMSB cases 
respectively. If the upper value of the top quark mass, M t = 182.3 GeV, is used, and an 
estimated 4 GeV theoretical uncertainty on Mh is included, one arrives at the conservative 
upper bound of Mh ^152 GeV for the value of the lighter Higgs boson in the unconstrained 
MSSM. Such a bound has important implications for searches and phenomenology of the 
MSSM Higgs bosons. 

Note added 

During the final stages of the preparation of the present work, a paper appeared [HE] , where 
a two-loop computation of the leading contributions to the MSSM Higgs boson self-energies 
for arbitrary external momentum is presented. In the numerical examples presented in 
Ref. [69J, the two-loop momentum-dependent corrections to the lighter Higgs boson mass 
Mh are of the order of a few hundred MeV, while the corresponding corrections to the heavier 
Higgs boson mass Mh are considerably smaller. These results are in good agreement with 
our empirical estimates given in section 13.3.31 We plan to include the results of Ref. (HO] in 
our spectrum calculators as soon as they become available in a suitable form. 
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Appendix A: EWSB conditions and Higgs boson masses 

In this appendix we provide for completeness the general formulae for the computation of 
the radiative corrections to the EWSB conditions and to the Higgs boson masses in the 
MSSM, following the notation of Ref. ^H]- We also explain how to adapt the formulae of 
Refs. |21I23 123 123 to the same notation. 

Once the requirement of correct EWSB is imposed, providing input values for the Z 
boson mass and tan/? = v 2 jv\ [where v% (z = 1, 2) are the vacuum expectation values of the 
neutral Higgs fields], the minimisation conditions on the MSSM effective potential translate 
into conditions for the parameters /i 2 and B: 
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where mz is the running mass for the Z boson defined in eq. (J2J), m|. (i = 1, 2) are the soft 
SUSY-breaking mass terms for the two Higgs doublets and U are the tadpole corrections, 
corresponding to diagrams with one incoming Higgs field. The squared physical masses for 
the CP-odd Higgs boson A and for the charged Higgs boson H ± can be computed as: 

Ml = 2B/sm2(3-ReU AA (Ml)+sm 2 (3— + cos 2 /3 — , (A3) 

v\ v 2 

M 2 H± = M 2 +M^ + Re [U AA (M 2 A )-U H+H -(M 2 H± ) + U^ W (M^)} , (A4) 

where M w is the physical W boson mass, and H AA (M\) , U H + H - (M^-±) and n^ w (M 2 F ) are 
the self-energies for the A , H ± and W bosons, respectively, each computed at a squared 
external momentum equal to the squared physical mass of the corresponding particle. 

Finally, the squared physical masses for the CP-even Higgs bosons, M 2 and M#, corre- 
spond to the real part of the poles of the 2x2 propagator matrix, i.e. the two solutions of 
the equation: 



Dot 



p 2 I-M 2 (p 2 ) =0, (A5) 



where the radiatively corrected mass matrix for the CP-even Higgs fields is: 
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In the above formula, m 2 A = 2 Bj sin 2(3 is the squared running mass for the A boson as 
defined in eq. (|Z|), and Uij(p 2 ) (i,j = 1,2) are the self-energies of the CP-even components 
of the Higgs fields. In Sof tSusy, SPheno and SuSpect an iterative procedure is employed 
for solving eq. (jA5|) . The Higgs mixing angle a, instead, is defined in the codes as the angle 
that diagonalises M 2 {p 2 ) for p 2 = M\. 

The explicit formulae for the one-loop parts of the self-energies and tadpole diagrams 
appearing in eqs. (jAlj) - ()A6j) can be found in Ref. ^Bj- In addition, SoftSusy, SPheno 
and SuSpect employ the formulae given in Refs. [2"31 12"51 l2"o1 12*5] for the leading two-loop 
corrections to the tadpoles and the neutral Higgs boson masses. 

The computations in Refs. [23 123 123 I2H] are performed in the effective potential ap- 
proach, which implies that the approximation of zero external momentum is used in the 
two-loop self-energies. Also, the formulae of Refs. [231231211121] assume that the tree-level 
mass matrix for the CP-even Higgs bosons is expressed in terms of the physical A-boson 
mass, Ma, whereas in eq. (|A6)1 the mass matrix is expressed in terms of the running mass 
m A . It is thus necessary to rearrange the formulae of Refs. (231231211121] for the corrections 
to the Higgs tadpoles, Sj, to the CP-odd Higgs mass, Am^, and to the CP-even Higgs mass 
matrix, (AM. 2 S ) , in order to adapt them to the notation of Ref. |lfij . In particular, the 
following equalities hold: 

U 

Si = -- 

Vi 



' (2 = 1,2), 



Am 2 A = -ReU AA (0) + sin 2 (3 — + cos 2 (3 — , 

vi v 2 



Am 2 A sin 2 /5+(A7W|)^ = -R e n n (0) + — , 
-Am 2 A smf3cos(3+(AM 2 s Y S = -Ren 12 (0) 



rr J. 

Am 2 A cos 2 p+ (AM 2 S Y = -Ren 22 (0) + — . 



Appendix B: Formulae for the (9(a&a T ) corrections 

We present in this appendix the explicit analytical formulae for the two-loop 0(ab<y T ) correc- 
tions to the MSSM Higgs masses and tadpoles in the DR renormalisation scheme. We follow 
very closely the lines of Refs. [23 123 123 123 121], to which we refer the reader for further 
details. The following results complete the computation of the two-loop corrections to the 
Higgs masses and tadpoles controlled by the third-family Yukawa couplings (the inclusion 
of the 0(a 2 ) corrections was discussed in Ref. [2E])- 

If we express the MSSM effective potential V e s in terms of DR-renormalised fields and 
parameters, the two-loop contribution to V c r involving both the bottom and tau Yukawa 
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couplings hb and h T reads, in units of iV c /(i67r 2 ) 



h b h T 

^ V br — y S 2 g b S 2 g T C(p b -(p T 



(Bl) 

where: c^ = cos0 and s^ = sin0; m? (q = T,b) are the field-dependent squark masses; 0g 
is the field-dependent squark mixing angle, defined in such a way that < 9 q < -k/2 (to be 
contrasted with the usual field-independent mixing angle 6q, such that — n/2 < 9q < vr/2); 
(p q is the phase in the off-diagonal element of the squark mass matrix. For the explicit 



Higgs field dependence of these parameters, see Refs. (221 E5j- Finally, the DR-subtracted 
two-loop integral J, with Q 2 being the renormalisation scale, is defined as 

Hmlml) = ml m 2 2 (l - In ^ (l - In ^ , (B2) 

Exploiting the field-dependence of the squark masses and mixing angles, we can express 
the tau/bottom corrections to the CP-even Higgs mass matrix, (AM|) C , to the Higgs 
tadpoles, Ej, and to the CP-odd Higgs mass, Am\, (see appendix A and Ref. J2E] for the 
various definitions) as 

ff 1 
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+ h T h b A T A b s 29T s 29b F 5 + 4:h T h b m T m b F 6 , (B3) 

-(AM|) 12 = h 2 T ^m T s 2 e T F 2 T + -h 2 T A T ^s 2 2dT Fl + h T h b m b iis 2dT Fl 
+ h 2 h \im b s 29b F 2 + -hlAbfi s 2 2f)b F 3 6 + h T h b m T /j, s 29b F b 
+ -h T hbS 2eT s 2 e b n(Ab + A r )F 5 , (B4) 

(AM 2 ^ = ]-h 2 T ^ 2 s 2 2eT Fl+ l -h 2 b ^ 2 s 2 2eb Fl + h T h bl ?s 29T s 2eb F^ (B5) 



v\lh. = m T A T s 29T F T + 2mlG T + m b A b s 2db F b + 2m 2 b G b , (B6) 

—v\ E 2 = m T /i tan (3 s 29t F t + m b // tan (3 s 29b F b , (B7) 

Am 2 , = -— (2h b h T F A - f^\ IT- k }^ Ab 2 F b \ (BS) 

c /3 Sj9 \ "^ - r< m^ - mf 2 
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In the equations above, A T and A\, are the soft supersymmetry-breaking trilinear couplings 
of the Higgs fields to the stau and sbottom squarks, /i is the Higgs mass term 11 in the 
superpotential, tan/3 = v 2 /vx is the ratio of the two Higgs vacuum expectation values and 
s 28 q (<? = r , b) refer to the usual field-independent squark mixing angles. The functions 
F? (i = 1,2,3,4) , F 5 , Fq , F q , G q and Fa are combinations of the derivatives of AVb T 
with respect to the field-dependent parameters, computed at the minimum of the effective 
potential; their explicit expressions are, in units of iV c /(i67r 2 ) 2 : 
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11 Note that we employ here the opposite convention for the sign of fi with respect to Ref. 26 , which 
explains the minus signs appearing in eqs. (|B4() , (|B7|I and 
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replacement r <-» 6. The function J g (g = r, b) entering the above equations is defined as: 



2 \ / 2 

vn%, \ r, i mi 



The results in Eqs. (|B9Jl - (|B17|) are valid when the one-loop 0(a T ) and 0(atb) corrections are 
written in terms of DR parameters computed at the renormalisation scale Q 2 . In alternative, 
an on-shell renormalisation scheme for the tau/stau and bottom/sbottom sectors that avoids 
the appearance of tan /^-enhanced terms in the two-loop part of the results should be devised 
along the lines of Refs. 
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